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The  objective  of  this  program  is  to  develop  and  validate  structural  de¬ 
sign  criteria  that  account  for  the  transient  overloads  produced  by  bird  and 
ice  impacts  on  turbine  engine  first-stage  f an/compressor  blades.  This  pro¬ 
gram  is  part  of  a  continuing  effort  of  the  Air  Force  to  minimize  ownership 
costs  by  placing  added  emphasis  on  the  derivation  of  durable,  damage-tolerant, 
advanced-turbomachinery  component  designs. 

Foreign  object  damage  (FOD)  in  aircraft  engines  has  been  an  increasing 
drain  on  defense  economic  resources.  Present-generation  engine  blading  de¬ 
signs  have  substantial  material  toughness  and  cross  section  to  meet  the 
mechanical  and  aeromechanical  requirements,  but  they  were  not  specifically 
designed  for  damage  tolerance.  The  prevailing  design  state  of  the  art  con¬ 
cerning  damage-tolerant  blading  is  mostly  empirical  and  based  on  experience 
from  major  incidents  and  associated  ad  hoc  testing.  The  empirical  approach 
to  designing  damage-tolerant  blading  requires  considerable  investment  of 
resources  and  is  sometimes  incapable  of  providing  directly  predictable  foreign- 
object-impact  response  of  blading.  This  approach  is  not  adequate  for  newer, 
more  damage-prone  blading  designs  incorporating  lightweight  materials  (in¬ 
cluding  advanced  composites)  and  thinner  cross  sections  to  achieve  improved 
performance . 


Comprehensive  design  criteria  based  upon  transient-stractur al-response 
tools  are  needed  to  achieve  more  efficient,  damage-tolerant  blading.  This 
program  will  provide  the  necessary  computer  tools  and  validation  testing  to 
establish  reliable  design  criteria  to  provide  direct  assessment  of  blade 
impact-damage  tolerance  and  identify  areas  for  improvement. 


This  program  consists  of  11  tasks  that  progressively  develop  the  FOD 
design  criteria  from  computer  models,  structural  element  and  material  pro¬ 
perty  tests,  and  static  and  rotating  single-blade  tests  to  full-  or  partial1 


stage  rotating 

tests.  The  11  tasks  are  identified  as  fol 

Task 

I 

- 

Design  System  Structure 

Task 

II 

- 

Transient  Response  Analysis  Models 

Task 

III 

- 

Impact  Loading  Model 

Task 

IV 

- 

Material  Response  and  Failure  Criteria 

Task 

V 

- 

Parametric  Analysis 

Task 

VI 

- 

Structural  Element  Tests 

Task 

VII 

- 

Error  Band  Analysis 

Task 

VIII 

- 

Foreign  Object  Impact  Design  Criteria 

Task 

IX 

- 

Single-Blade  Impact  Tests 

Task 

X 

- 

Full-Stage  Response  Prediction 

Task 

XI 

- 

Full-Scale  Impact  Tests 

The  design  analysis  methods  and  failure  criteria  derived  in  the  course 
of  this  program  will  be  applicable  both  to  advanced  composite  materials  and 
to  monolithic  materials  for  current  and  advanced  f an/compressor  blading. 
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% 
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2 . 0  SUMMARY 


The  program  consists  of  11  tasks  to  establish  reliable  FOD  design  crite¬ 
ria.  The  emphasis  will  be  on  determining  the  effects  of  the  ingestion  of 
birds  and  ice  during  takeoff  and  landing  (which  present  the  greatest  hazards 
to  aircraft  engines). 

The  11  tasks  progressively  develop  the  FOD  design  criteria  from  transient 
structural-response  computer  models,  test  specimens,  and  single-blade  tests  to 
full-  or  partial-stage  rotating  rig  tests.  The  transient-response  computer 
models  will  be  based  on  existing  mathematical  approaches  and  on  existing  com¬ 
puting  programs  that  will  be  adjusted  as  required  for  the  specific  application 
of  this  study. 

Task  I  has  been  completed.  A  design  structure  was  developed  for  the 
purposes  of  providing  a  consistent  framework  in  which  to  constitute  the  design 
criteria,  guiding  the  overall  technical  effort,  and  formulating  a  data  founda¬ 
tion  to  define  the  problem  environment. 

Tasks  II  through  IX  are  currently  underway  and  address  the  development 
of  the  transient-response  and  loading  models,  local  and  structural  impact 
testing,  material-property  determination,  error-band  analysis,  and  the  formu¬ 
lation  of  design  criteria. 

The  preliminary  and  final  designs  of  the  transient-structural-response 
computer  programs  have  been  developed,  and  checkout  is  nearing  completion. 

In  addition,  the  coding  for  interfacing  the  bird-loading  model  (developed  by 
the  Univeristy  of  Dayton  Research  Institute)  with  the  final-design  response 
model  (developed  by  General  Electric)  has  been  completed  and  checked  out. 

Models  of  the  three  first-stage  blades  (representing  airfoils  from  Air 
Force  inventory,  production-development,  and  advanced-development  engines) 
have  been  constructed,  and  bench  impact  testing  of  these  blades  has  been 
completed . 

Composite  and  metallic  specimens  for  local  impact  and  structural-element 
bench  impact  and  material-property  tests  have  been  designed  and  fabricated. 
Plans  defining  the  testing  of  these  specimens  and  real  blades  have  been 
completed  and  approved  by  the  Air  Force.  Real-blade  impact  testing  has  been 
completed  to  establish  the  maximum  strain  rates  needed  for  the  material  pro¬ 
perty  tests. 

A  parametric  matrix  has  been  formulated  to  define  the  conditions  and 
geometries  that  will  be  analyzed  in  Tasks  V  and  X  and  to  provide  guidance 
in  the  structural  tests  of  Tasks  VI  and  IX.  Transient  analyses  have  been 
performance  on  olate-like  structural  models. 
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Work  was  performed  on  the  compilation  of  relevant  blade  and  impactor 
properties  for  use  in  the  error-band  analysis.  FOD  design  criteria  were  es¬ 
tablished  within  the  framework  of  the  projected  capabilities  of  the  transient 
structural-analysis  computer  tools  developed  in  Tasks  II  and  III. 


Single-blade  impact  testing  (whirligig)  was  completed  on  the  J79  com- 
oressor  Stage  1  steel  blade,  and  data  reduct  ion/analysis  is  underway.  This 
blade  represents  Air  Force  inventory  airfoils. 


3.0  TASK  I  -  DESIGN  SYSTEM  STRUCTURE 


A  design  system  structure  was  developed,  through  the  Task  I  effort,  to 
establish  a  consistent  framework  in  which  to  constitute  FOD  design  criteria 
and  guide  the  overall  technical  effort.  This  design  system  structure  defines 
the  objective  and  interfaces  of  the  various  multidisciplined  tasks  that  must 
be  satisfied  to  achieved  meaningful  design  criteria,  and  it  includes  a  data 
foundation  that  describes  the  problem  environment. 

Figure  1  depicts  the  interrelations  between  the  major  distinct  activi¬ 
ties  of  quantification  of  experience,  analytic  development,  bench  tests, 
rotating  rig  tests,  and  the  work  flow  of  the  program. 

The  results  of  the  Task  I  effort  were  documented  in  a  Research/Develop¬ 
ment  Test  Plan^  that  was  approved  by  the  Air  Force. 

To  conduct  the  program  described  in  Figure  1,  a  team  of  technical  per¬ 
sonnel  has  been  gathered  from  several  organizations  to  utilize  their  specific 
technical  expertise  and  FOD  experience.  These  organizations  are: 

•  General  Electric  Aircraft  Engine  Group  (AEG) 

•  General  Electric  Corporate  Research  and  Development  Center  (CR&DC) 

•  University  of  Dayton  Research  Institute  (UDRI) 

Figures  2  and  3  show  the  functional  responsibilities  of  the  participants  and  the 
organization . 


■^Superscript  numbers  refer  to  references. 
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•  Coordinate  Program 

•  Parametric  Analyses 

»  Rotating  Single-Blade 
and  Full-Stage  Testin 


4.0  TASK  II  -  TRANSIENT  RESPONSE  ANALYSIS  MODELS 


The  first-  and  second- level  model s^  selected  for  use  in  the  transient- 
response  analysis  of  first-stage  fan  and  compressor  blades  were  approved  by 
the  Air  Force. 


4.1  INTRODUCTION 


The  transient  response  due  to  foreign  object  impact  on  first-stage  fan 
and  compressor  blades  will  be  analyzed  on  two  different  levels  during  the 
course  of  this  program.  The  rationale  for  this  approach  is  based  on  the 
fact  that  each  level  of  analysis  possesses  intrinsic  advantages  which  make 
it  very  attractive  for  the  study  of  different  aspects  of  the  FOD  problem. 

The  first  level  of  analysis  will  be  accomplished  with  a  modified  version 
of  the  finite-element  program  NONSAP.  This  analysis  will  be  capable  of 
carrying  out  a  very  detailed  analysis  embracing  local  impact  damage  as  well  as 
gross  structural  damage.  NDNSAF  is  a  very  general  and  accurate  package  for 
nonlinear,  finite-element  analyses  and  employs  three-dimensional  (3-D)  elements. 
The  modified  version  used  in  this  program  will  be  called  NQSAPM.  This  pro¬ 
gram  has  a  library  of  three  different  time- integration  techniques  (Wilson 
Newmark,  central  difference)  as  well  as  two  elastoplast ic  hardening  models 
(kinematic  hardening  and  isotropic  hardening).  This  versatility  will  prove 
useful  in  finding  the  most  efficient  strategy  for  solution  of  the  dynamic 
problem.  However,  the  detail  and  versatility  of  this  NONSAP  package  is 
obtained  at  the  expense  of  high  computing  costs  which  could  prove  to  be  a 
problem  during  parametric  design  studies.  As  a  result,  the  second-level 
analysis  based  on  S.  Levy's  component-element  method  of  transient  dynamics^ 
will  be  employed  in  situations  where  cost  becomes  an  important  factor.  The 
savings  in  the  Component  Element  Method  (COMET)  are  accomplished  through  the 
use  of  simplified  element  models  and  an  explicit  time- integration  technique 
that  allows  nonlinear  structural-dynamics  problems  to  be  solved  without 
physically  assembling  stiffness  matrices.  This  Level  2  analysis  will  be 
somewhat  less  detailed,  but  experience  in  other  fields  including  pipe  whip\ 
locomotive  dynamics,  and  nuclear  fuel  rod  vibration^  shows  that  the  technique 
is  an  accurate  as  well  as  efficient  tool.  In  addition  to  being  of  use  in 
parametric  studies,  the  lower  cost  of  this  level  of  analysis  should  prove  to 
be  an  advantage  in  blade-interaction  studies. 

No  dynamic  crack-growth  analyses  will  be  incorporated  into  the  first- 
and  second-level  analyses.  Stresses  and  strains  will  simply  be  calculated 
and  then  used  to  indicate  damage  based  on  the  design  criteria  to  be  formed  as 
part  of  Task  VIII.  Plastic  curl-back  of  blades  will  be  considered.  Graphical 
output  of  stresses  and  strains  will  be  provided,  and  a  method  of  implementing 
the  results  of  the  impact-damage  design  criteria  within  the  structural  codes 
will  be  incorporated.  Based  on  the  transient-response  models  of  Task  II  and 
the  damage  criteria  formulated  in  Task  VIII,  the  final  package  would  then 
indicate  probable  areas  of  failure  as  they  are  encountered  in  the  analysis. 
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4.2  FIRST-LEVEL  RESPONSE  MODEL  (NOSAPM) 


As  indicated  above,  the  program  NOSAPM  contains  a  wide  range  of  tools 
that  will  be  necessary  in  modeling  the  FOD  problem.  Specifically,  in  this 
analysis,  the  3-D  isoparametric  element  will  be  used.  This  element  can  be 
used  with  a  variable  number  of  nodal  points  ranging  from  8  to  21  (see  Figure 
4).  It  is  anticipated  that  16-noded  elements  will  be  used  to  allow  the  blades 
to  be  modeled  with  one  element  through  the  thickness  while  still  including  the 
bending  effects  which  will  be  most  important  in  this  problem.  NOSAPM  is 
capable  of  large-displacement,  elastoplastic  analysis  and  uses  an  efficient, 
out-of-core,  equation  solver.  The  inelastic  material  behavior  will  be 
modeled  with  a  bilinear  stress-strain  curve,  and  both  isotropic  and  kinematic 
hardening  will  be  available  for  application.  In  the  case  of  composite  mate¬ 
rials,  a  linear  orthotropic  material  model  is  already  available.  Centrifugal 
stiffening  has  been  incorporated  in  NOSAPM  through  the  addition  of  a  body- 
force  (centrifugal)  load  vector.  This  load  vector  is  used  to  calculate  cen¬ 
trifugal  displacements  which  are  used  as  initial  conditions  for  the  transient 
response.  The  constant  centrifugal-load  vector  is  also  added  to  the  impact 
loads  for  each  time  step.  The  centrifugal  displacements  calculated  by  the 
initial  state  of  the  blade  before  the  impact  are  also  utilized  in  the  calcu¬ 
lation  of  eigenfrequencies  and  eigenmodes  that  indicate  the  effects  of 
centrifugal  stiffening. 

If  the  strain-rate  effects  are  found  to  be  substantial,  in  the  course  of 
material  characterization,  they  will  also  be  included.  Flexible  root  fixity 
as  well  as  blade  interaction  at  the  shrouds  will  be  handled  with  spring  ele¬ 
ments.  Table  1  is  a  summary  of  the  capabilities  of  NOSAPM. 


Table  1.  First-Level  (NOSAPM)  Analysis  Based  on  Finite-Element  Method. 


•  Blade  model  formulated  from  3-D  solid  elements  of  from  8  to  2 1 
nodes . 

•  Sixteen-noded  elements  have  been  selected  (midsided  nodes  only 
for  in-plane  directions)  for  the  blade  models. 

•  NOSAPM  is  capable  of  large-displacement,  elastoplastic  analysis 
and  uses  an  efficient,  out-of-core,  equation  solver. 

•  The  inelastic  material  behavior  is  modeled  with  a  bilinear 
stress-strain  curve,  and  both  isotropic  and  kinematic  hardening 
are  available  for  application. 

•  For  composite  materials,  a  linear  orthotropic  material  model  is 
used . 

•  Centrifugal  stiffening  has  been  incorporated  through  the 
addition  of  a  body-force  (centrifugal)  load  vector. 

•  Flexible  root  fixity  as  well  as  blade  interaction  at  the  shrouds 
is  handled  with  spring  (Truss)  elements. 
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Figures  5,  6,  and  7  show  NOSAPM  three-dimensional  models  of  the  three 
first-stage  blades  selected  for  analysis  in  Task  V.  Figure  5  shows  the  un¬ 
shrouded,  stainless  steel,  J79  compressor  first-stage  blade.  Figure  6  shows 
the  unshrouded,  boron/ aluminum  (B/Al)  composite,  first-stage  fan  blade  de¬ 
veloped  in  the  Advanced  Propulsion  System  Integration  Program  (APSI).  And 
Figure  7  shows  the  shrouded,  titanium,  F101  first-stage  fan  blade. 

In  order  to  facilitate  the  running  of  NOSAPM  from  Evendale,  a  data-deck- 
generator  program  has  been  written.  This  program  generates  card  images  of 
the  NOSAPM  input  data  based  on  element  and  node  files  (created  with  an  exis¬ 
ting  mesh-generator  program  available  at  Evendale)  and  on  user  input  specified 
in  the  interactive  time-sharing  mode.  The  input  data  decks  so  created  are 
transmitted  to  the  WPAFB  6600  systems  via  the  remote  Evendale  input/output 
(I/O)  site  to  establish  data  files.  Control  of  job  execution  in  the  batch 
mode  is  then  accomplished  at  Evendale,  through  time-sharing,  with  the  job 
output  directed  to  the  remote  I/O  site  at  Evendale.  Work  is  proceeding  on 
the  coding  of  a  deck-generator  program  for  the  second-level  analysis,  and 
Figure  8  shows  a  schematic  of  the  steps  required  in  generating  and  running 
blade  models  with  the  first-  and  second-level  analyses. 

Figure  9  shows  the  results  of  a  check  case  used  to  verify  the  accuracy 
of  the  NOSAPM  solution  for  displacements  induced  by  centrifugal  forces. 

The  figure  indicates  excellent  agreement  between  a  closed-form  solution  and 
NOSAPM  for  centrifugal  displacements  in^the  radial  direction  for  a  flat- 
plate  rotating  at  5000  rpm.  In  NOSAPM  this  flat  plate  was  modeled  by  the 
thirty  16-noded-element  configuration  shown  in  Figure  10.  Referring  to 
Figure  11,  the  expression  for  the  closed-form  solution  was  developed  as 
follows.  The  elemental  centrifugal  force  is: 

dF  =  A  dX  (p/g)  w2  (ri  +  X)  (1) 

where  p/g  is  the  mass  density  and  w  is  che  angular  velocity.  The  centri¬ 
fugal  stress  developed  by  the  plate  mass  lying  above  the  section  defined  by 


Z  is: 

o  -  F/A  =  (p/g)  w2  zJ£(ri  +  X)  dX  (2) 

a  =  (p/g)  i2  [r^  ( l  -  Z)  +  ( i2  -  Z2)/2]  (3) 

Radial  strain  is: 

ez  =  dW/dZ  (4) 

Radial  displacement  is: 

w=rzezdz  (5) 

J0 

w  =  ( pu>2  /gE)  o/Zf-ia  -  Z)  +  (£2  -  Z2)  /  2  )  dz  (6) 

W  =  (pw2/gE)  [ri  UZ  -  Z2/2)  +  U2Z  -  Z3/3)/2]  (7) 
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Figure  6.  Three-Dimensional,  Finite-Element 
Model  of  APSI  B/Al  Fan  Blade. 
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Figure  8.  Blade  Modeling  and  Analysis  System. 
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Figure  9.  Radial  Displacement  Due  to  Centifugal 


Significant  progress  has  been  made  in  the  integration  of  plotting 
capabilities  into  the  first-level  NOSAPM  analysis  to  aid  data  interpretation. 
The  graphics  package  MOVIE-BYU  is  being  applied  to  represent  displacement  and 
strain  data  from  NOSAPM.  MOVIE-BYU  is  a  3-D,  finite-element,  plotting  code 
in  general  use  today;  it  offers  a  variety  of  options.  Preprocessing  of  data 
output  by  NOSAPM  should  provide  all  the  information  required  for  this  program. 
Briefly,  these  are  the  capabilities  currently  available  with  MOVIE-P.yU  and 
the  associated  preprocessing  completed  to  date: 

1.  Deformed  mesh  plots  can  be  obtained  at  predetermined  time 
interval s . 

a.  These  plots  can  be  made  with  the  option  of  suppressing 
hidden  lines  -  thus  improving  the  ease  of  visualization. 

Figures  12  through  1A  illustrate  typical  results  from  a 
normal  impact. 

b.  The  distortion  of  the  mesh  can  be  scaled  to  any  size  the 
user  desires  in  order  to  make  small  displacements  visible. 
Figure  15  illustrates  this  option.  Here,  the  displacements 
in  Figure  13  are  amplified  fivefold. 

c.  The  user  can  move  his  field  of  view  of  the  structure  and 
"zoom  in"  to  study  local  behavior. 

2.  Plots  of  the  effective  strain  on  the  surface  of  the  elements  can 
be  displayed  as  contours. 

a.  These  plots  can  be  displayed  either  on  a  deformed  or  an 
undeformed  structure. 

b.  With  proper  supporting  equipment,  these  contours  can  be 
translated  into  color  plots.  This  makes  it  possible,  in 
turn,  to  identify  surface  areas  which  have  surpassed  yield 
strain  by  coding  these  areas  all  to  be  a  particular  color  - 
red  for  example. 

3.  A  method  of  obtaining  surface  strains  on  the  deformed  structure 
has  been  formulated  for  the  preprocessor  and  is  being  programmed 
now. 


A. 3  SECOND-LEVEL  RESPONSE  MODEL  (COMET-BLADE) 


Work  on  the  development  of  the  second-level  response  model  is  basically 
complete.  Essentially,  the  same  physical  effects  are  considered  in  the  Level 
2  (COMET-BLADE)  analysis  as  in  the  Level  l,  but  the  element  approach  is  some¬ 
what  simplified  to  gain  efficiency.  The  Level  2  analysis  is  based  on  the 
component-element  method^  and  provides  an  analytical  computer  capability  for 
predicting  gross  and  local  structural  response  of  blades  to  representative 
large- bod y-impac t  loading  stimuli. 
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Bending  and  membrane  effects  ace  modeled  with  a  system  of  plate  elements. 
These  elements  couple  membrane  and  bending  forces  and  are  similar  to  those 
that  have  ben  highly  successful  in  modeling  both  plasticity  and  hysteresis  in 
the  General  Electric  COMET-PIPE  ccxnputer  program^. 

The  formulation  and  programming  of  the  plate  elements  for  large  deforma¬ 
tions,  material  yielding,  and  hysteresis  effects  are  complete.  These  elements 
have  nine  nodes  and  address  three  cases:  (a)  an  element  around  a  central  node 
for  interior  regions  of  the  blade,  (b)  a  corner  element,  and  (c)  a  middle-of- 
an-edge  element.  The  COMET  plate  element  has  nine  nodes  and  is  shown  in  Figure 
16;  Figure  17  illustrates  Level  2  modeling. 

Level  2  analysis  also  includes  centrifugal  loading  effects,  in  a  manner 
similar  to  that  of  the  first-level  analysis,  and  provides  for  the  specifica¬ 
tion  of  damping  that  can  be  proportional  to  mass,  to  stiffness,  or  to  a  com¬ 
bination  of  the  two.  The  blade  can  have  twist  and  camber  and  can  vary  in 
thickness.  Boundary  conditions,  such  as  restraint  of  rotation  and  displace¬ 
ment,  can  be  specified.  An  elastic,  anisotropic  option  is  provided,  and 
spring-damper  elements  can  be  used  to  restrain  nodes  on  the  blade. 

The  blade  is  modeled  by  giving  the  positions  of  nodes  in  global  coordi¬ 
nates,  the  blade  thickness,  the  material  real-stress  versus  real-strain  curve, 
the  density,  the  damping,  the  speed  in  rpm,  and  the  boundary  conditions.  The 
loads  can  be  applied  at  any  node  in  any  direction.  Step-by-step  numerical 
integration  is  used  to  obtain  a  solution  for  the  response. 

Table  2  is  a  summary  of  COMET  capabilities. 


4.4  INTERFACE  WITH  THE  LOADING  MODEL 


The  interface  between  the  response  and  loading  models  will  be  described 
in  detail  in  the  loading  model  section  (Section  5).  A  brief  description  of 
the  interface  follows. 

Since  loading  due  to  a  foreign- obj ect  impact  will  depend  upon  the  defor¬ 
mation  of  the  blade,  interaction  between  the  structural  model  and  the  impactor 
model  will  be  necessary.  This  requires  the  interfacing  of  the  response 
routines  developed  by  CR&D  with  the  loading  model  developed  by  UDRI. 

The  interface  between  the  first-level  response  model  developed  by  CR&D 
and  the  loading  model  developed  by  UDRI  has  been  accomplished.  Geometric  and 
dynamic  information  is  periodically  transferred  from  the  structural  analysis 
and  used  to  update  the  impact  loads.  This  need  not  be  done  at  every  struc¬ 
tural  time  step  of  the  transient  analysis  since  large  time  steps  may  be  suit¬ 
able  without  loss  of  accuracy. 

The  interface  with  a  loading  model  developed  by  CR&D  has  also  been  in¬ 
cluded  in  the  second-level  response  model.  Figure  18  shows  the  linkup  between 
the  loading  model  subroutine  and  the  executive  or  main  response-model  program. 
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Table  2.  Second-Level  (COMET-BLADE)  Analysis 
Based  on  Component-Element  Method. 


The  component- element  method  provides  an  efficient,  nonlinear, 
analysis  capability  that  does  not  require  the  assembly  of  a 
global  stiffness  matrix  or  the  solution  of  simultaneous  equations 

Simplified  element  models  and  an  explicit  time-integration  tech¬ 
nique  are  combined  with  an  approach  that  treats  external  (bird 
loading)  and  internal  restoring  forces  as  forces  externally 
applied  to  nodal  points  representing  a  distributed  mass  system. 

No  iteration  is  necessary,  and  the  solution  for  future  displace¬ 
ments  and  strains  is  obtained  by  calculating  the  current  values 
of  the  nodal  accelerations  using  the  current  values  of  all  ex¬ 
ternally  applied  forces. 

Z^i  =  Qi/mi  =  current  acceleration  at  node  point  i 

ZLi  =  2Z®i  -  Z~^i  +  Z^i  (At)^  =  future  displacement 

Centrifugal  stiffening  is  incorporated  through  a  body-force  load 
vec  tor . 

The  blade  model  is  constructed  in  3-D  space  from  a  system  of  nine 
noded  plate  elements  and  addresses  twist,  camber,  and  thickness 
variation. 

The  plate  elements  model  bending  and  membrane  forces  and  include 
the  effects  of  large  displacements,  material  yielding,  and 
hysteresis . 
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Figure  18.  Linkup  Between  the  Loading  Model  Subroutine  and 
the  Executive  or  Main  Response-Model  Program. 


5.0  TASK  III  -  IMPACT  LOADING  MODEL 


5.1  INTRODUCTION 


The  objective  of  this  task  is  to  generate  an  impact-loading  model,  for 
birds  and  ice,  which  can  be  interfaced  and  used  with  the  transient-response, 
structural-analysis  tools  developed  in  Task  II.  The  effort  in  this  task  is 
divided  into  four  subtasks:  (a)  simulated  bird  and  ice  impact  parameters, 
(b)  loading-model  sensitivity  study,  (c)  formulat ion/ interf ac ing  of  loading 
model  with  structural  models,  and  (d)  experimental  verification  of  loading 
model . 


5.2  SUBTASK  A  -  SIMULATED  BIRD  AND  ICE  IMPACT  PARAMETERS 

The  work  on  this  subtask  has  been  completed.  A  report^  written  by  UDRI 
documenting  the  effort  has  been  completed  and  transmitted  to  the  Air  Force. 
This  report  presents  the  results  obtained  from  spatially  and  temporally  re¬ 
solved  pressure  measurements  made  during  impacts  of  artificial  (gelatin)  birds 
and  ice  on  rigid  targets.  The  data  indicate  that  gelatin  behaves  the  same  as 
real  birds  during  impact  and  that  the  analytically  predicted  pressure  distri¬ 
butions  obtained  with  the  loading  model  of  Subtask  C  agree  with  the  measured 
pressures  both  for  bird  impacts  and  for  ice  impacts. 

5.3  SUBTASK  B  -  LOADING  MODEL  SENSITIVITY  STUDY 

The  work  under  this  subtask  has  been  completed,  and  a  UDRI  report  on  the 
effort  is  being  reviewed  prepratory  to  transmittal  to  the  Air  Force.  The 
purpose  of  the  sensitivity  study  is  to  provide  specific  guidance  as  to  which 
aspects  of  the  loading-model  formulation  must  be  modeled  most  accurately  to 
properly  predict  specific  aspects  of  blade  response. 

The  Materially  And  Geometrically  Nonlinear  Analysis  (MAGNA)  finite- 
element  code  developed  by  UDRI  was  used  for  the  numerical  analysis  in  t in¬ 
sensitivity  study.  The  complex  geometry  of  a  real  fan  blade  was  avoided  in 
this  study  by  using  a  flat,  cantilever  plate.  This  simplification  of  real 
bird-impact  conditions  in  a  jet  engine  permitted  a  less  ambiguous  study 
of  the  loading  details.  A  flat  plate  7.62  cm  (3  in.)  across  the  chord  by 
24. 46  cm  (9.63  in.)  span  by  3.81  cm  (0.15  in.)  thick  was  chosen  as  repre¬ 
sentative  of  medium-sized  jet  engine  fan  blades.  The  geometry  is  shown  in 
Figure  19.  The  material  properties  of  Type  410  stainless  steel  were  em¬ 
ployed  in  the  model. 

The  finite-element  model  of  the  cantilever  plate  is  shown  in  Figure  20. 
Three-dimensional  solid-brick,  isoparametric  elements  with  21  nodes  (the  21st 
node  was  located  at  the  center  of  the  loaded  face)  were  used  to  model  tin¬ 
plate.  The  additional  node  on  the  loaded  face  permitted  more  accurate  load 
definition  over  the  loaded  elements.  Analyses  were  performed  to  account  for 
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geometric  and  material  nonlinearities.  The  nonlinear  dynamic  analysis  re¬ 
quires  considerable  computing  time  for  small  time  steps  where  a  long  response 
calculation  is  desired.  The  time  step  for  most  analyses  was  fifteen  micro¬ 
seconds.  This  short  time  step  dictated  by  solution  accuracy  made  only  early 
time,  local  response  analysis  practicable  for  the  large  number  of  analyses 
performed.  The  local  response  (as  opposed  to  overall,  late-time  deflections) 
was  determined  to  be  more  valuable  as  a  measure  of  sensitivity  to  loading 
function  details.  The  overall,  maximum,  late-time  deflection  of  the  canti¬ 
lever  plate  was  always  directly  proportional  to  the  projectile  momentum  at 
impact.  Thus,  either  bird  size  or  impact  velocity,  not  the  details  of  spa¬ 
tial  or  temporal  load  distribution,  affected  the  overall  deflection. 

5.4  SUBTASK  C  -  FORMULATION/ INTERFACING  OF  LOADING  MODELS  WITH 

STRUCTURAL  MODELS 

5.4.1  Introduct ion 

The  loads  exerted  on  a  structure  during  an  impact  can  be  characterized 
in  terms  of  a  number  of  fundamental  parameters.  The  most  important  of  these 
parameters  are  the  total  impulse  imparted  to  the  structure,  the  forces  (i.e., 
the  rate  of  impulse  transfer),  the  pressures  (i.e.,  the  spatial  distribution 
of  the  force),  and  the  time  variation  of  forces  and  pressures. 

Loading  models  have  been  developed  that  address  the  parameters  identi¬ 
fied  above  to  the  extent  required  to  accurately  define  bird  and  ice  loading 
on  a  target  which  not  only  translates  and  rotates  during  impact  but  also  de¬ 
forms  locally.  The  loading  models  compute  the  pressure  distribution  exerted 
on  surfaces  during  bird  and  ice  impacts  (soft-body  impacts).  The  bird  and 
ice  impacts  are  considered  as  fluid  dynamic  in  nature  and  are  modeled  as 
fluid  jets  impinging  on  arbitrarily  shaped,  three-dimensional,  deformable 
surfaces . 

A  quasi-steady,  potential-flow  analysis  is  applied  to  the  jet-impact 
problem  thereby  reducing  it  to  the  problem  of  solving  Laplace's  equation. 

The  surface-singularity  technique  is  used  to  solve  Laplace's  equation. 

A  subroutine  for  computing  pressure  distributions  on  either  rotating  or 
nonrotating  turbine  engine  components  has  been  developed.  A  model  for  treat¬ 
ing  slicing  of  birds  and  ice  by  rotating  blades  is  incorporated  into  this 
routine.  The  loading-model  routine  is  specifically  developed  to  be  inter¬ 
faced  with  finite-element,  structural-analysis,  computer  programs.  Given  the 
instantaneous  impacted  surface  shape  and  displacement  velocity,  the  loading- 
model  program  computes  the  pressure  distribution  existing  on  the  impacted 
surface.  The  shape  of  the  impact  surface  is  fed  to  the  loading  model  in  the 
form  of  finite-element,  surface-nodal-point  coordinates.  From  this  informa¬ 
tion,  the  loading-model  computer  program  constructs  a  system  of  quadrilateral 
elements  covering  the  impacted  surface.  A  uniform  distribution  of  fluid  dy¬ 
namic  sources  of  unknown  but  determinable  strength  is  applied  to  each  of  these 
quadrilateral  elements.  The  strength  of  each  source  distribution  is  deter¬ 
mined  from  boundary  condition  considerations.  The  jet-velocity  field  and  the 


pressure  distribution  are  ultimately  determined  from  the  velocity  field  in¬ 
duced  by  the  distribution  of  sources  over  the  impact  surface  combined  with 
the  jet-velocitv  field  which  would  exist  without  the  presence  of  the  impacted 
sur  face . 


5.4.2  The  Slicing  Model 

Because  the  forces  required  to  slice  a  soft  body  are  small  compared  to 
the  forces  required  to  decelerate  the  body,  the  slicing  forces  can  be  ignored. 
Therefore,  the  slicing-model  development  reduces  to  simply  a  geometric  prob¬ 
lem  of  determining  the  dimensions  and  weight  of  a  slice.  For  birds  and  ice 
spheres  (such  as  hail),  the  velocity  of  the  ingested  object  is  small  compared 
to  the  velocity  of  the  aircraft;  thus,  the  object  velocity  can  be  ignored. 

For  slab  ice,  such  as  ice  breaking  loose  from  an  engine  nacelle,  the  velocity 
of  the  slab  relative  to  the  nacelle  is  not  well-defined,  but  we  can  expect 
this  velocity  to  be  small  in  comparison  to  the  aircraft  velocity. 


5.4. 2.1  Slicing  Model  Development  for  Birds 


The  bird  is  idealized  as  a  right  circular  cylinder  with  a  length-to- 
diameter  ratio  of  2,  and  the  velocity  of  the  bird  relative  to  the  aircraft 
is  taken  to  be  equal  and  opposite  to  the  aircraft  velocity.  In  the  following 
analysis,  a  coordinate  system  attached  to  the  blade  is  used.  The  bird/blade 
interaction  geometry  in  such  a  coordinate  system  is  shown  in  Figure  21.  The 
following  information  is  assumed  to  be  known  (supplied  as  input  data  to  the 
loading  model  subroutine): 


N  -  number  of  blades  per  stage 
n  -  blade  rotational  speed  (rpm) 


-  distance  from  the  rotational  axis  of  the  rotor  to  the  point  on 
the  blade  at  which  the  center  of  the  impact  occurs 

6  -  blade  orientation  angle  (defined  as  the  angle  6  in  Figure  21) 


Vjj  -  axial  velocity  of  the  bird 

Wfj  -  bird  weight 

Pb  -  bird  density  (0.91  g/cm^  or  0.033047  lbm/in.^ 

An  infinite  number  of  slice  shapes  are  possible  for  a  given  set  of  these 
input  parameters;  the  shape  depends  on  the  orientation  of  the  bird  relative 
to  the  blade  and  on  the  span  location  on  the  blade  at  which  the  impact  oc¬ 
curs.  Since  worst-case  slice  shapes  are  desired  (i.e.,  slice  shapes  having 
the  largest  possible  slice  mass),  the  orientation  of  the  axis  of  the  bird 
and  the  center  of  impact  of  a  slice  are  chosen  to  produce  a  slice  having  the 
largest  possible  mass.  The  largest  slice  mass  occurs  when  the  axis  of  the 
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slice  is  coincident  with  the  axis  of  the  right  circular  cylinder.  The  slice 
geometry  depicted  in  Figure  21  corresponds  to  this  worst-case  situation. 

With  the  bird  idealized  as  a  right  circular  cylinder  with  L  *  2D,  the 
diameter  of  the  bird  is  determined  from  the  bird  weight  and  density  as: 

D  =  (2Wb/irpb)  !/3  (8) 

The  tangential  velocity  of  the  bird,  Vt,  is  computed  at  the  impact 
radius  Z£  and  is  given  by: 

Vt  =  2irnZi/60  (9) 

And  the  magnitude  of  the  velocity  of  the  bird  relative  to  the  blade,  Vr ,  is 
given  by: 


Vr  =  /(Vb2  +  Vt2)  (10) 

The  direction  of  Vr  relative  to  the  plane  of  rotation  of  the  blade  at  the 
impact  location,  Z^,  is: 

0  =  sin"1  (Vb/Vr)  (11) 


Vr  is  assumed  to  have  no  component  in  the  radial  direction.  The  impact 
angle,  defined  as  0  in  Figure  21,  is: 

0=6-0  (12) 

The  blade  spacing  (defined  at  Z^  and  denoted  by  S  in  Figure  21)  is  given  by: 

S  =  2wZi/N  (13) 


The  bird  slice  width,  h,  is  then  given  by: 
h  =  S  sin  0  =  ( 2wZ £/N)  (Vb/Vr) 

Referring  to  Figure  21,  we  find  that  the  bird-slice 
+h/2 


Wsb  Pb^ 


-h/2 


we  ight 


is  given  by: 


( 14) 


(15) 


wsb  =  2pbD  [h/(D/2)2  -  (h/2)2  +  (D2/2)  sin"1  (h/D )] 


(16) 


35 


The  use  of  Equations  9  and  10  in  Equation  14  shows  that  the  bird-slice 
width  at  any  impact  radius  depends  only  on  the  blade  parameters  (N  and  n)  and 
the  aircraft  speed.  The  bird-slice  weight  at  any  impact  radius  depends,  in 
addition  to  the  above  three  parameters,  only  on  the  density  and  total  weight 
of  the  bird. 

If  h/sin  0  is  greater  than  the  chord  of  the  blade  at  Z^,  then  not  all 
of  the  slice  will  impact  on  the  blade. 


5. 4. 2. 2  Slicing  Model  for  Ice  Spheres 


For  ice  spheres  (such  as  hail),  the  axial  velocity  of  the  sphere  rela¬ 
tive  to  the  blade  is  taken  to  be  equal  and  opposite  to  the  aircraft  velocity, 
just  as  for  birds.  The  expressions  for  slice-width,  impingement  angle,  and 
Vr  are  the  same  as  for  birds.  The  slice  weight,  Ws^,  given  by: 


h/2 


Wsis  "  2lroi  f  KD/2)2  ~  52] 


dC 


or: 


Ws.s  =  2*Pi  [(D/2)2  (h/2)  -  (h/2)3/5] 


(17) 


(18) 


where  h/2  £  D/2,  and  pj^  is  the  ice  density.  For  ice  spheres  it  is  presumed 
that  the  diameter  of  the  sphere  is  a  specified  quantity  along  with  the  ice 
density. 


5.4. 2. 3  Slicing  Model  for  Ice  Slab 

For  slabs  of  ice,  the  operator-specified  information  is  presumed  to  be 
the  length  of  the  slab  ( L ) ,  the  thickwess  of  the  slab  (&Z)  (i.e.,  how  much  of 
the  blade  span  is  to  be  exposed  to  the  slab),  and  the  axial  velocity  of  the 
slab  relative  to  the  nacelle  'this  replaces  Vj,  in  Equations  10,  11,  and  14). 
The  slice  width,  h,  is  again  given  by  Equation  14  with  V^,  replaced  by  the 
axial  velocity  of  the  slab  relative  to  the  nacelle.  The  slice  mass  is  then 
given  by 

wss  =  Pi<h)(AZ)L  (19) 

5.4.3  Bird-Slice  Geometry  Parametric  Study  Results 

Using  the  equations  developed  in  Section  5.4.2. 1,  the  bird-slice  parame¬ 
ters  were  computed  both  for  starling  impacts  (3-ounce  birds)  and  for  big-bird 
impacts  (1.5-pound  birds).  Three  different  blade  configurations  (J79,  APSI, 
and  F 1 0 1 )  were  considered,  and  typical  values  of  rotor  speed,  number  of  blades 
per  stage,  blade  orientation  angle,  and  impact  radius  were  used.  The  input 
parameters  and  the  computed  quantities  are  shown  for  the  two  cases  of  bird 
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weights  in  Tables  3  and  4.  The  primary  observation  from  this  parametric  study 
is  that  bird  (and  ice  sphere)  impacts  are  highly  oblique. 


5.4.4  Fluid  Flow  Model  of  Impact 

It  has  been  demonstrated  in  Reference  6  that  for  birds  having  an  L/D  of 
2  or  greater  the  largest  contributing  factor  in  the  total  impulse  of  a  bird 
impacting  a  target  is  the  steady-flow  phase  of  the  impact.  Thus,  in  the  cur¬ 
rent  version  of  the  loading  model,  the  pressure  distribution  is  based  on  a 
steady-flow  analysis.  This  method  of  analysis  is  extended  into  a  dynamic 
analysis,  for  deforming  targets,  by  treating  the  flow  at  any  instant  of  time 
as  a  quasi-steady  flow  using  the  instantaneous  relative  velocity  between  the 
bird  slice  and  the  deforming  target  as  the  characteristic  velocity.  The  in¬ 
stantaneous  shape  of  the  target  is  used  to  define  the  surface  on  which  impact 
occurs . 

In  the  analysis  of  the  impact  process,  the  strength  of  the  projectile 
material  is  considered  to  be  negligible.  This  assumption  is  quite  reasonable 
for  the  typical  projectile  and  target  materials  of  interest. 

We  assume  that,  at  least  as  a  first  approximation,  the  fluid  flow  may  be 
treated  as  incompressible.  This  is  a  reasonable  assumption  because  the  mea¬ 
sured  steady-state  pressures  are  quite  small  in  comparison  to  the  pressures 
required  to  produce  significant  density  changes  in  birds  or  ice. 

By  virtue  of  the  incompressible  flow  and  steady  flow  assumptions  and  the 
negligible  strength  of  the  projectile  material,  the  problem  of  predicting 
soft-body  impact  loads  is  amenable  to  analysis  within  the  following  real  con¬ 
straints.  The  computing  time  required  to  obtain  a  reasonable  solution  must 
be  small  in  comparison  to  the  computing  time  required  to  compute  the  struc¬ 
tural  response,  and  the  computer  storage  requirement  for  the  loading  model 
must  be  small  in  comparison  to  the  structural-response-analysis  computer  pro¬ 
gram  storage  requirements. 

Three-dimensional  potential  flow  theory  was  chosen  as  the  most  appropri¬ 
ate  method  for  modeling  the  impact  process.  Our  initial  work  in  using  this 
approach  was  quite  successful^.  In  that  work,  we  used  the  method  of  surface 
singularities  to  determine  the  velocity  and  pressure  fields  caused  by  circular 
jets  impacting  at  oblique  angles  on  flat,  rigid  plates.  This  same  technique 
was  used  to  develop  the  loading  model.  The  method  developed  earlier  was  gen¬ 
eralized  to  include  the  impact-surface  curvature  and  deformation  velocity, 
an  arbitrary  cross-sectional  area  of  the  impacting  flow,  and  a  method  for 
generating  planar  quadrilateral  surface  elements  from  given  finite-element, 
surface-nodal-point  locations. 
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Table  4.  Blade  Calculation  Input  Data  Sheet  for  1.5-lb  Bird 
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5 . 4 . 4 .  1  Surface  Singularity  Method  -  Superposition  of  Onset  Flow 
and  Surface  Source  Distribution 

The  impacting  bird  slice  is  idealized  as  an  onset  flow,  of  velocity 
Vr,  uniformly  distributed  over  the  cross-sectional  area  of  the  slice.  The 
fluid  velocity  at  a  point  is  expressed  as  the  negative  gradient  of  a  potential 
function  that:  satisfies  Laplace's  equation  in  the  region  exterior  to  the 
blade,  has  a  zero  normal  derivative  on  the  blade  surface,  and  approaches  the 
proper  uniform  stream  potential  function.  The  flow  is  viewed  as  an  onset 
flow  potential  and  disturbance  potential  due  to  the  blade.  The  disturbance 
potential  satisfies  Laplace's  equation  in  the  exterior  and  the  appropriate 
normal  derivative  condition  on  the  surface  and  vanishes  at  infinity.  From 
potential  theory,  it  can  be  shown  that  this  potential  may  be  evaluated  in 
terms  of  a  surface  source  density  distribution  with  which  the  blade  surface 
may  be  considered  to  be  covered.  The  incorporation  of  the  normal  derivative 
condition  on  the  surface  produces  a  source  density  distribution,  over  the 
surface,  described  by  a  two-dimensional  Fredholm  integral  equation  of  the 
second  kind.  Once  this  equation  is  solved  for  the  source  densitv  distribu¬ 
tion,  first  the  disturbance  potential  and  then  the  disturbance  flow  veloci¬ 
ties  may  be  obtained. 


5.4.4. 2  The  Body  Surface  Approximation 

The  body  surface  must  be  suitably  approximated  in  the  numerical  solution 
of  the  Fredholm  integral  equation  for  the  surface  source  density  distribution. 
The  body  surface  is  approximated  by  a  large  number  of  smaii,  plane  elements 
formed  from  the  original  points  defining  the  body  surface  The  calculations 
employ  plane,  quadrilateral  elements;  the  source  density  is  assumed  constant 
over  each  of  these  elements.  With  this  assumption,  the  problem  of  determining 
the  continuous  function  for  the  source  density  reduces  to  one  of  determining 
a  finite  number  of  values  -  one  for  each  of  the  planar  elements. 


5.4.4. 3  Solution  Procedure 

This  section  presents  a  brief  description  of  the  method  of  solution. 

As  a  first  step,  it  is  necessary  to  input  a  set  of  points,  in  three- 
dimensional  space,  describing  the  body  surface.  This  input  should  be  chosen 
to  provide  the  best  representation  of  the  body  with  the  fewest  possible 
points.  Four  of  these  inputs  are  employed  to  generate  one  plane,  quadri¬ 
lateral,  source  element  of  the  body  surface. 

Once  the  values  of  the  surface  source  density  have  been  determined,  the 
velocity  components  at  each  null  point  are  calculated,  and  the  appropriate 
components  of  the  onset  flow  are  added.  The  resultant  velocities  and  pres¬ 
sures  are  easily  computed  from  the  velocity  components. 
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5.4.5  Structure  of  the  Loading  Model  and  Interface  With  The 
First-Level  Response  Model 


We  have  incorporated  in  our  analysis  a  procedure  for  finite-difference 
numerical  computation  of  the  potential  flow  about  arbitrary,  three-dimensional 
bodies;  this  procedure  has  certain  features  that  render  it  compatible  with  the 
first-level,  finite-element,  structural-response-analysis  code.  Specifically, 
the  nodal  points  in  the  finite-element  program  become  approximate  locations  of 
the  null  point  (defined  as  the  point  where  the  quadrilateral  element  induces 
no  velocity  in  its  own  plane)  around  which  the  plane,  quadrilateral  elements 
are  constructed  in  the  potential-flow  analysis.  This  greatly  facilitates  the 
coupling  of  the  finite-difference  loading  model  and  the  finite-element  struc¬ 
tural  code. 

In  loading  used  for  any  arbitrarily  shaped  impact  surface,  the  impact 
area  is  divided  into  small,  flat  elements;  a  uniform  distribution  of  sources 
is  assumed  to  cov£r  each  area.  At  the  beginning  of  impact,  an  initial  pres¬ 
sure  distribution  is  computed  for  the  undeformed  blade.  During  an  impact  in 
which  local  deformation  takes  place,  the  deformed  shape  of  the  impact  zone 
is  calculated  in  the  dynamic  structural-response  analysis.  After  signifi¬ 
cant  deformation  has  occurred,  the  geometry  of  the  impact  zone  is  provided  to 
the  loading  model.  The  loading  model  is  then  used  to  calculate  a  new  pres¬ 
sure  distribution.  As  the  structural-analysis  calculation  proceeds,  the 
local  shape,  the  location,  and  the  velocity  of  the  impact  area  are  updated 
and  passed  to  the  loading  model  at  appropriate  time  intervals.  The  loading 
model,  in  turn,  provides  updated  pressure-distribution  information  for  the 
structural-response  computation.  The  loading  model  is  fully  interactive  with 
the  structural-response  calculation.  The  duration  of  the  impact  is  computed 
by  keeping  track  of  how  much  slice  has  been  consumed  during  each  time  inter¬ 
val  . 


A  summary  of  the  interface  between  the  loading  and  response  models  fol¬ 
lows  : 

•  The  response  model  provides  displacements  and  velocities  at  the 
surface  nodes  of  the  structural  elements  in  the  impact  region. 

These  structural  elements  are  then  subdivided  into  smaller  elements 
corresponding  to  the  finer  grid  of  the  loading  model  as  specified 
by  user  input,  and  the  structural-element  shape  functions  are  then 
used  to  compute  the  displacements  and  velocities  at  the  nodes  of 
the  smaller  elements.  These  displacements  and  velocities  are  pro¬ 
vided  to  the  loading  model  at  specified  time  steps. 

•  The  slice-size  calculations  are  based  on  user  input  and  are  per¬ 
formed  by  the  response  program.  The  impactors  available  to  the 
user  include  the  following: 

Bird  of  any  weight  represented  by  a  right  circular  cylinder 
having  L/D  ratio  equal  to  2.0. 

Ice  ball  of  any  diameter. 
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Slab  ice  of  any  thickness  and  length. 


User-provided  information  (in  addition  to  the  geometry  and  material 
properties  of  the  finite-element  response  model)  will  consist  of 
the  velocity  of  the  aircraft,  the  rotational  velocity  of  the  fan  or 
compressor,  the  number  of  blades  in  the  stage,  the  spanwise  loca¬ 
tion  of  the  impact,  and  the  specification  of  the  impactor. 

•  The  loading  model  provides  the  pressures  at  the  nodes  of  the  smalle 
elements,  and  the  structural-element  shape  functions  are  then  used 
to  calculate  the  pressures  at  the  nodal  points  of  the  structural 
elements.  These  pressures  are  then  converted  to  equivalent  nodal 
forces  through  a  consistent  force  formation. 

•  A  flag  is  returned  by  the  loads  model  when  the  impactor  has  been 
consumed  (i.e.,  when  it  has  traveled  its  own  length).  When  this 
occurs  the  response  model  will  no  longer  call  the  loading  model 
subroutine,  and  the  nodal  pressures  are  set  to  zero. 

The  loading  model  is  capable  of  detailed  interaction  with  the  structural 
response  model  and  of  dealing  with  target  translation,  rotation,  and  local 
deformation.  The  load/response  coupling  modeled  in  this  formulation  should 
be  capable  of  accurately  predicting  both  overall  target  response  and  local 
de  format  ion . 

The  principal  limitation  of  the  loading  model  is  that  it  does  not  in¬ 
clude  transient  effects.  The  most  significant  of  these  is  the  shock  effect; 
however,  the  porosity  present  in  birds  appreciably  reduces  the  shock  stresses 
without  significantly  affecting  steady  Elow  pressures.  In  addition,  impact 
obliquity  reduces  the  relative  importance  of  shock  stresses.  Therefore,  it 
is  not  obvious  that  neglect  of  the  shock  aspects  of  bird  impact  on  blades  is 
a  significant  deficiency.  Another  transient  aspect  of  bird  impacts  is  the 
variation  of  slice  cross  section  at  the  target  surface  during  impact.  This 
variation  results  mainly  in  a  time  variation  of  the  impact  area.  If  the  flow 
remains  quasi-steady  during  these  variations  (i.e.,  the  "velocity"  of  the 
variation  is  low  with  respect  to  local  sound  speed  -  probably  a  good  assump¬ 
tion),  then  the  loading  model  can  be  modified  to  describe  these  effects.  The 
size  and  geometry  of  the  "jet"  that  flows  onto  the  target  surface  must  be  up¬ 
dated  incrementally  to  describe  the  variation  with  time  of  the  impact  area. 

5.4.6  Structure  of  the  Loading  Model  and  Interface  With  The 
Second-Level  Response  Model 

A  simplified,  interactive,  load  model  was  developed  by  CR&D  for  the 
second-level  response  code.  This  model  is  less  elaborate  than  the  potential- 
flow  approach  incorporated  by  UDRI  and  is  based  in  part  upon  experimental 
impact  information  available  in  Reference  6.  However,  it  does  offer  several 
advantages;  as  a  result,  it  has  been  incorporatsd  in  the  Level  2  stiuctural- 
response  model.  Among  these  advantages,  perhaps  the  most  important  is  sim¬ 
plicity,  The  pressure  calculations  in  this  model  require  only  the  geometry 
and  the  velocity  at  the  point  for  which  pressure  is  to  be  calculated;  as  a 
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result,  no  linear-equation  solution  is  required.  This  reduces  the  time  for 
calculation.  In  addition,  the  experimental  basis  for  the  model  allows  for 
the  inclusion  ot  the  initial  pressure  shock  during  impact. 

The  general  operation  of  this  load  model  is  quite  simple.  The  bird  and 
the  impact  geometry  are  exactly  the  same  as  those  defined  and  applied  in 
UDRI's  model.  The  pressures  are  calculated  at  the  structural-node  points. 

Each  time  the  loads  are  calculated,  the  program  determines  which  nodes  of  the 
structure  are  within  the  impact  zone.  It  also  identifies  the  location  of  the 
trailing  end  of  the  impacting  bird  with  respect  to  the  current  blade  location. 
Pressure  loading  is  terminated  when  the  end  of  the  bird-slice  passes  the  de¬ 
formed  location  of  the  blade.  During  the  impact  event,  the  pressures  at  nodes 
which  fall  within  the  impact  zone  are  calculated  as  follows. 

With  knowledge  of  the  structural  geometry  and  nodal  velocities,  as  well 
as  the  bird-slice  velocity,  a  relative  velocity  between  bird  and  structure 
normal  to  the  plane  of  the  blade  is  defined.  This  velocity  is  then  used  to 
calculate  an  average  pressure  based  upon  the  simple  impulse-momentum  con¬ 
siderations  shown  in  Figure  22.  This  pressure  reflects  a  dependence  upon  the 
structural  deformation  occurring  during  impact. 

Use  is  now  made  of  two  experimental  observations.  The  first  observation 
is  made  with  respect  to  Figure  14  in  Reference  6.  This  figure  reflects  the 
observation  that  the  loading  on  the  blade  rises  to  a  peak  of  twice  the  average 
pressure  at  approximately  20%  of  total  impact  time  before  decreasing  to  zero. 
Figure  23  shows  the  temporal  variation  of  pressure  that  is  nond imens ional ized 
with  respect  to  the  average,  steady-state,  dynamic  pressure.  These  experi¬ 
mental  data  are  incorporated  in  the  loading  model  as  a  time-scaling  factor: 
FCTRT . 

The  second  observation  is  made  with  respect  to  Figure  34  of  Reference  6. 
For  oblique  impacts,  this  information  indicates  that  the  pressure  is  maximum 
at  the  leading  edge  and  approximately  equal  to  twice  the  average,  steady- 
state,  dynamic  pressure.  The  pressure  decreases  to  zero  at  a  given  distance 
from  the  leading  edge.  Refer  to  Figure  24  for  a  description  of  the  spatial 
variation  of  pressure.  This  behavior  is  incorporated  in  the  model  through 
the  use  of  a  second  scaling  factor,  FCTRG,  which  is  also  applied  to  the  aver¬ 
age  dynamic  pressure.  Thus,  the  pressure  applied  at  any  given  node  of  the 
structure  is: 

P  =  ( FCTRG) ( FCTRT ) ( p V2 / g )  sin2  9  (20) 

The  calculation  is  simple,  interactive,  and  representative  of  the  experimental 
results  reported  in  Reference  6. 

5.4.7  Summary 

Loading  models  have  been  developed  for  soft-body  impacts  on  structures 
that  respond  to  the  impact  in  a  very  general  manner.  The  loading  models  were 
developed  from  a  fluid-dynamic  description  of  a  soft  body  impacting  (flowing 
onto)  a  structure.  The  structure  can  translate,  rotate,  and  deform  locally 
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Impulse.  I  =  FAt  =  MAV  =  Momentum  Transfer 
F  =  I/At  =  MAV /At 
F  =  MV  sin  9  (V/£) 

F  =  p  A  sin  9 

p  =  (.:/g)V^  sin^  9  =  Average  Steady-State  Dynamic  Pressure 
o  =  Bird  Density  =  0.033047  lbm/in.^ 

V  sin  9  =  Normal  Velocity  of  Bird  Relative  to  Blade 


Figure  22.  Impulse-Momentum  Principle;  Bird  is  Treated  as  a  Fluid  Body. 


during  the  impact.  The  loading  model  for  the  first-level  response  code  is 
based  on  an  approximate  solution  to  the  problem  of  the  quasi-steadv  flow  of  a 
jet  of  arbitrary  cross  section  flowing  onto  a  surface  with  an  arbitrary  shape. 

The  cross-sectional  area  of  the  jet  is  determined  from  an  analysis  of  the 
slicing  of  soft  bodies  by  a  row  of  fan  or  compressor  blades.  The  loading- 
model  pressure  distributions  for  the  first-level  response  code  show  excellent 
agreement  with  pressure  distributions  obtained  in  experiments  in  which  birds, 
substitute  birds,  and  ice  spheres  impacted  on  rigid  targets.  The  loading 
model  has  been  successfully  interfaced  with  the  first-level,  finite-element, 
structural-analysis  computer  program. 

The  overall  validity  and  limitations  of  the  loading  models  will  be  deter¬ 
mined  by  comparing  the  structural  response  predicted  by  the  combined  loading 
model / struc tural-ana lys is  model  with  that  observed  in  experiments. 


5.5  SUBTASK  D  -  EXPERIMENTAL  VERIFICATION  OF  LOADING  MODELS 


In  order  to  check  the  adequacy  of  the  integrated  loading  and  response 
models,  comparisons  will  be  made  between  the  analysis  predictions  and  the 
test  results  for  impacts  on  cantilevered  plates.  A  range  of  impact  parame¬ 
ters  (mass  for  substitute  birds  and  ice,  velocity,  and  incidence)  will  be 
employed  to  ensure  a  full  range  of  specimen  response  modes  (i.e.,  elastic  and 
elastic-plastic  structural  deformation  and  severe  local  deformation).  The 
predictions  of  the  complete  model  will  be  compared  to  experimentally  derived 
results  obtained  from  Task  VI  -  Structural  Element  tests. 

The  experimental  work  has  been  completed.  Displacement  data  collected 
with  the  Moire  fringe  device  for  cantilever  plates  struck  by  simulated 
(gelatin)  birds  has  been  reduced.  These  measurements  provide  detailed,  out- 
of-plane-displacement  data  at  several  increments  of  time  and  are  presented 
in  Figure  25.  The  excellent  agreement  between  these  measurements  and  analysis 
predictions  is  discussed  in  Section  7.4.3  of  this  report. 
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Figure  25.  Deflected-Shape  Plots  of  Cantilever  Plate  Struck  by  Simulated  (Gelatin)  Bird  for 
Anal yt i cal -Model  Validation. 


6.0  TASK  IV  -  MATERIAL  RESPONSE  AND  FAILURE  CRITERIA 


6.  I  SUBTASK  A  -  GROSS  STRUCTURAL  DAMAGE  PROPERTIES 

6.1.1  Introduction 

Testing  has  been  performed  in  two  areas  to  establish  gross  structural- 
material  response  and  failure  data.  First,  bench  impact  tests  have  been 
conducted  using  the  full-scale  blades  selected  in  Task  V.  The  strain  rates 
pertinent  to  the  problem  of  foreign-ob j ec t  impacts  have  been  determined  from 
these  tests.  Impact  tests,  using  the  parameters  of  Tables  3  and  4,  have 
been  conducted  on  three  different  type  blades,  two  of  which  are  unshrouded. 

The  two  unshrouded  blades  are  the  J79  compressor  first-stage  stainless  steel 
blade  and  the  APSI  first-stage  fan  B/Al  composite  blade.  These  blades  were 
cantilever-mounted  at  the  root.  The  third  blade  type  is  the  F101  first-stage 
fan  titanium  blade.  This  blade  is  tip  shrouded  and,  therefore,  was  fastened 
at  both  ends.  Each  blade  type  was  impacted  at  the  30  and  70%  span  locations, 
and  two  blades  of  each  blade  type  were  used. 

In  the  second  part  of  Subtask  A,  tests  are  being  conducted  to  obtain 
material-property  data  for  the  relevant  strain  rates.  Forged  bar  stock 
was  used  to  make  the  metal  specimens  (8A1-1MO-1V  titanium  to  represent  the 
F10 1  blade  and  410  stainless  steel  to  represent  the  J79  blade).  The  410  is 
being  substituted  for  403  stainless  steel  which  is  unavailable.  Chemically, 
the  two  alloys  are  the  same  except  for  a  small  amount  of  chromium;  this 
element  does  not  affect  the  room- temperature  mechanical  properties.  Standard, 
cylindrical  specimens  are  used  for  the  tension  tests  at  low  and  intermediate 
strain  rates,  and  the  entire  stress-strain  curve  is  obtained  using  strain 
gages  and  an  extensomer.  The  change  in  diameter  is  used  to  calculate  Poisson's 
ratio.  The  high- strain-rate  tests  are  conducted  on  the  split  Hopkinson  bar 
using  the  standard  tension-test  specimen  for  that  apparatus. 

The  mechanical-property  tests  have  been  completed  for  the  metallic  mate¬ 
rials,  and  testing  is  in  progress  to  establish  the  mechanical  properties  for 
the  composite  material. 

The  advanced  composite  selected  for  the  study  is  boron/ aluminum ,  used  in 
the  APSI  Stage  1  fan  blade,  with  a  fiber  layup  of  (0/22/0/-22 ) .  The  4-mil 
boron  fibers  are  in  a  matrix  of  2024  aluminum.  Each  ply  is  4.7  mils  thick 
and  does  not,  in  general,  run  the  length  of  the  blade.  The  center  ply  of  the 
blade  is  a  plate  of  aluminum.  On  either  side  of  it  is  a  ply  of  stainless 
steel  wire  mesh.  The  next  plies  on  both  sides  are  the  B/Al  layup.  The  entire 
outer  blade  surface  is  covered  with  a  ply  of  stainless  steel  wire  mesh.  The 
leading  edge  has  a  thin  coat  of  pure  nickel. 

The  material  properties  are  needed  in  three  principal  directions  because 
the  composite  is  anisotropic.  The  data  needed  are  the  tensile  stress-strain 
curves  for  each  direction,  three  Poisson's  ratios,  and  three  shear  stress- 
strain  curves.  However,  the  collection  of  this  information  is  simplified  con¬ 
siderably  by  the  use  of  laminate  theory.  This  theory  employs  the  results  of 
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material-property  tests  on  unidirectional  specimens  to  calculate  the  bulk 
properties  of  the  angle- ply  layup. 


6.1.2  Real-Blade  Strain  Rates 

Testing  of  the  full-scale  blades  to  obtain  the  maximum  strain  rates  has 
been  completed.  The  highest  strain  rates  developed  were  less  than  400  in./ 
in. /sec  in  any  of  the  types  of  blades  tested  (J79,  F101,  and  APSI).  For  the 
J79  impacts,  the  highest  strain  rate  for  the  30%  span  impact  was  344  in./in.- 
sec  for  a  1.5-pound-bird  impact.  This  value  was  for  a  gage  located  on  the 
midchord  at  the  root  in  the  spanwise  direction.  For  the  70%  span  impact  of 
the  J79  blade,  the  highest  strain  rate  was  362  in. /in. /sec  for  a  1.5-pound- 
bird  impact.  This  value  was  for  a  gage  directly  behind  the  impact  site  and 
oriented  in  the  chordwi.se  direction. 

The  highest  strain  rate  for  the  30%  span  impacts  of  the  tip- re strained 
F101  blade  was  341  in. /in. /sec.  For  the  70%  span  impacts  of  this  blade,  the 
highest  strain  rate  was  367  in. /in. /sec.  Both  strain  rates  were  for  1.5-pound- 
bird  impacts  and  for  a  gage  located  at  the  blade  tip,  along  the  leading  edge, 
in  the  spanwise  direction. 

For  the  APSI  blade,  the  highest  strain  rate  developed  was  366  in./ in. /sec 
for  a  70%  span  impact  of  a  3-ounce  bird.  This  occurred  at  a  gage  at  the  mid¬ 
chord  of  the  root  in  the  spanwise  direction.  The  highest  strain  rate  for  the 
30%  span  impact  was  362  in. /in. /sec  for  a  3-ounce  bird.  This  gage  location 
was  behind  the  impact  site,  and  the  orientation  of  the  gate  was  in  the  direc¬ 
tion  of  the  chord. 


6.1.3  Mater ial-Propert5 


The  material-characterization  tests  on  the  two  metals,  410  stainless 
steel  (410  substituted  for  403)  used  for  the  J79  compressor  Stage  1  blade  and 
titanium  8AI-IM0-IV  used  in  the  F101  fan  Stage  1  blade,  have  been  completed. 
The  tests  on  the  advanced  composite  (B/Al)  used  in  the  APSI  fan  Stage  1  blade 
have  not  yet  started.  As  mentioned  previously,  the  advanced  composite  used  in 
the  APSI  blade  is  actually  a  hybrid  composite.  This  blade  has  the  following 
layup:  SSWM/0°/22,’/0*/-22o/SSWM/Al/SSWM/-22'/0'’/22V0*/SSWM.  It  will,  there¬ 
fore,  be  necessary  to  characterize  the  stainless  steel  wire  mesh  also. 


Tensile  tests  conducted  at  the  strain  rates  shown  in  Table  5  provide  the 
desired  mechanc ial- property  data  for  the  metallic  materials.  This  range  of 
rates  encompasses  the  range  of  strain  rates  measured  at  select  locations  on 
he  full-scale  blade  tests,  12-370  str ain/second . 
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Table  5.  Strain  Rates  of  Metallic  Material 
Characterization  Tests. 


Material 

Strain  Rates  ( strain/ seconds) 

8AI-IM0-IV  Titanium 

403  Stainless  Steel 

10-3,  10"1,  1,  500 

10-3,  10_1,  1,  700 

The  low  and  intermediate  strain-rate  tests  (0.001  -  1.0  c/sec)  have  been 
conducted  per  ASTM  E8-69^),  and  the  specimen  definition  is  shown  in  Figure 
26.  A  standard  tensile  specimen  used  for  the  high-strain-rate  tests  on  the 
split  Hopkinson  bar  is  shown  in  Figure  27.  Table  6  shows  the  measured  densities 
of  9A1-1MO-1V  titanium  and  410  stainless  steel.  The  mechanical  properties 
acquired  from  the  various  strain-rate  tests  are  presented  in  Tables  7  and  8. 
Figures  28  and  29  are  plots  of  Young's  modulus  and  yield  stress  versus  strain 
rate . 


Table  6 .  Density. 


Material 

Density 

8AI-IM0-IV  Titanium 

410  Stainless  Steel 

4.32  g/ cm3 

7.57  g / cm3 

Relative  to  the  high-strain-rate  properties,  it  will  be  noted  that  a 
very  small  specimen  was  used  to  satisfy  assumptions  made  in  Hopkinson  bar 
theory  equations.  This  presents  problems  in  obtaining  accurate  data  from  the 
tests,  particularly  in  the  elastic  region.  One  problem  is  that  achievable 
machining  tolerances  on  the  specimen  in  relation  to  load-surface  alignment 
make  it  difficult  to  resolve  small  strains  accurately.  Secondly,  an  assump¬ 
tion  made  in  developing  the  Hopkinson  bar  equations  is  that  many  stress-wave 
reflections  occur  in  the  specimen  before  a  state  of  dynamic  equilibrium  is 
reached.  This  does  not  happen  until  some  plastic  deformation  has  taken 
place.  In  addition,  once  necking  occurs  in  the  specimen,  a  uniform  stress 
field  no  long  exists;  thus,  the  equations  for  stress  and  strain  are  invalid. 
This  is  a  significant  factor  for  the  small  specimens. 

However,  the  split  Hopkinson  bar  apparatus  represents  the  state-of-the- 
art  technique  for  obtaining  str ess- str ain  information  at  high  strain  rates. 
The  stress-strain  curves  are  considered  most  accurate  in  the  range  of  2  to 
10%  strain.  The  yield-stress  values  listed  in  Tables  7  and  8  represent  the 
point  on  the  curves  where  the  load  ceases  to  increase  significantly  for  a 
rather  large  increase  in  displacement.  (The  0,2%  offset  method  obviously 
cannot  be  used  because  of  the  lack  of  data  in  the  elastic  region.)  The  ulti¬ 
mate  stress  and  ultimate  strain  values  from  these  tests  were  the  final  point 
of  the  curve.  Although  the  values  from  the  high-strain-rate  tests  cannot  be 
used  as  exact  data  points,  they  are  useful  in  establishing  the  trend  of  the 
strain  rate  dependence  of  the  subject  materials  over  nearly  seven  decades  of 
strain  rate. 
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Threaded  Tension  Hopkinson  Bar  Specimen 


Diameter 

Diameter 

Diametei 


Figure  27.  Specimen  for  High-Strain-Rate  Tests. 


Table  7.  Mechanical  Properties  Obtained  for  8AI-IM0-IV  Titanium. 


Spec imen 

Strain 
Rate 
( e / sec) 

Elast ic 
Modulus 
(x  106  psi) 

Shear 
Modulus 
(x  10^  psi) 

Yield 

Stress 

(ksi) 

U1 1 imate 
Stress 
(ksi) 

Ultimate 

Strain 

(percent) 

1-Ti 

0.001 

16.9 

6.5 

148.9 

202.5 

28.2 

2-Ti 

0.001 

17.3 

6.7 

148.7 

200.8 

27.2 

3-Ti 

0.001 

17.4 

6.8 

148.7 

199.3 

21.8 

4-Ti 

0. 10 

15.8 

6.3 

156.3 

188.5 

13.7 

5-Ti 

0. 10 

16.9 

6.8 

153.0 

189. 1 

15.9 

6-Ti 

0. 10 

17.1 

6.8 

156.4 

194.5 

15.3 

7-Ti 

1.0 

17.3 

6.8 

164.0 

184.6 

15.7 

Ti-3 

550 

* 

* 

203.0 

161.1 

13.7 

Ti-2 

560 

* 

* 

202.0 

154.2 

14.5 

Ti-1 

580 

* 

* 

208.0 

158.7 

15.3 

Poisson's  Ratio  =  0.279 


^Elastic  properties  are  not  obtainable  from  the  high-strain-rate  tests. 
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Table  8.  Mechanical  Properties  Obtained  for  410  Stainless  Steel. 


Spec  imen 

St  r  a  in 
Rate 
( e/ sec ) 

Elast ic 
Modulus 
(x  10 ^  psi) 

Shear 

Mod  ul  us 
(x  10^  psi) 

Yield 

Stress 

(ksi) 

Ul  t imate 
Stress 
(ksi) 

Ul t imate 
Strain 
(percent) 

4-SS 

0.001 

27.4 

10.6 

84.2 

207.2 

26.7 

5-SS 

0.001 

30.  5 

11.7 

84.8 

207.2 

27.2 

6-SS 

0.001 

33.7 

12.9 

84.9 

211.9 

26.9 

9-SS 

0.  10 

27.7 

11.2 

87 . 2 

191.9 

25.7 

10-SS 

0.  10 

26.4 

10.2 

87.0 

188.3 

26.8 

11-SS 

0.  10 

28.2 

10.9 

89.2 

198.3 

23.5 

12-SS 

1.0 

27.2 

10.4 

86.6 

182. 7 

23.8 

13-SS 

1.0 

26.1 

10.0 

91.3 

185.5 

25.4 

14-SS 

1.0 

27.2 

10.4 

94.5 

183.8 

24.6 

15 -SS 

1.0 

28.0 

88.8 

184.7 

25.6 

SS-5 

680.0 

* 

* 

124.0 

129.5 

16.9 

SS-4 

760.0 

* 

* 

12  3.0 

130.0 

17.2 

SS-3 

700.0 

* 

* 

131.0 

13  1 . 4 

18 .0 

SS-2 

530.0 

* 

★ 

128.0 

irk 

~kk 

SS- 1 

340.0 

k 

* 

126.0 

irk 

— 

kk 

Poisson's  Ratio  =  0.297 


* 


Elastic  properties  are  not  obtainable  from  the  high-strain-rate  tests. 


These  specimens  did  not  fail. 
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Figure  28.  Young's  Modulus  Vs.  Strain  Rate 


ure  29.  Yield  Stress  Vs.  Strain  Rate 


The  test  matrix  that  will  be  used  for  the  characterization  of  the  boron/ 
aluminum  and  the  stainless  steel  wire  mesh  is  shown  in  Table  9.  For  the  ten¬ 
sile  tests  in  the  direction  of  the  symmetry  axes  (0°  and  90°),  a  modified 
ITTRI  specimen  configuration  will  be  used®.  This  specimen  is  shown  in  Figure 
30.  For  the  shear  tests  in  the  interlaminar  plane,  the  10°  off-axis  shear 
test  developed  by  Chamis  and  Sinclair^  will  be  used  for  the  B/Al.  A  draw¬ 
ing  of  the  specimen  appears  in  Figure  31.  Torsion  rod  specimens  will  be  used 
to  determine  the  transverse  shear  properties  for  both  the  B/Al  and  the  stain¬ 
less  steel  wire  mesh.  The  in-plane  shear  properties  for  the  stainless  steel 
wire  mesh  will  be  determined  with  the  ±45°  shear  panel  specimen. 

Table  9.  Characterization  Tests  for  B/Al  and 
Stainless  Steel  Wire  Mesh. 


Material  Test  Description  Strain  Rates 


B/Al  0°  Tension  1,10 

90°  Tension  1,10 

10°  Off-axis  Shear  10 

Rod  Torsion  Shear  10 

SSWM  Tension  (in  plane)  1,10 

Torsion  (through  the  thickness)  1,10 

±45°  Shear  10 

Rod  Torsion  Shear  10 


(  sec~ 1) 


The  through- the- thickness  tensile  specimen  shown  in  Figure  32  will  be 
used  to  determine  the  normal- to- the- plane  tensile  properties  for  the  stain¬ 
less  steel  wire  mesh.  The  tensile  specimens  (including  10°  off-axis  and  ±45° 
panels)  will  be  instrumented  with  high-resistance  strain  gages.  For  the 
through-the-thickness  and  torsion  rod  specimens,  it  is  planned  to  use  an  ex- 
tensometer  to  measure  displacement. 

It  is  not  anticipated  that  the  composite  will  be  significantly  strain- 
rate  sensitive.  Strain-rate  effects  in  B/Al  were  investigated  by  Krinke, 
Barber,  and  Nicholas^®  using  a  three-point  bend  (Charpy  impact)  technique. 
They  found  no  significant  strain-rate  effects  over  six  orders  of  magnitude; 
the  highest  strain  rate  employed  in  the  study  was  100  sec-*-. 


6.2  SUBTASK  B  -  LOCAL  LEADING-EDGE  DAMAGE 

6.2.1  Introduc  tion 

This  subtask  is  concerned  with  the  generation  of  empirical  design  data 
for  lopal  leading-edge  damage  to  first-stage  fan  and  compressor  blades. 
Most^df  the  testing  to  determine  and  quantify  the  damage  caused  by  leading- 
edge  impacts  for  a  range  of  pertinent  impact  conditions  has  been  accom¬ 
plished.  The  results  of  this  study  will  include  data  expressed  in  terms  of 
the  residual  properties  of  the  materials  which  are  directly  relevant  to  fan 
and  compressor  designers. 
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6.2.2  Overall  Approach 


Three  types  of  blade  materials,  geometries,  and  sizes  have  been  investi¬ 
gated  using  substitute  bird  materials,  ice,  granite  pebbles,  and  steel  spheres 
as  the  impactors.  This  range  of  impactors  is  representative  of  the  types  of 
objects  commonly  ingested  into  engines.  The  three  blade  types  are  the  J79 
blade  using  type  403  stainless  steel,  the  F101  blade  using  8-1-1  titanium, 
and  the  APSI  metal-matrix  boron/ aluminum  blade.  Fabricated  metal  and  compos¬ 
ite  specimens  and  actual  APSI  blades  have  been  tested.  The  geometries  of  the 
leading  edges  of  the  fabricated  specimens  are  similar  to  the  geometries  of 
the  blade  types  at  the  50%  span  location.  Three  specific  investigations  have 
been  conducted  in  this  subtask..  The  first  investigation  established  a  base  of 
damage  data.  The  second  addressed  scaling  and  subscale  tests.  The  third  in¬ 
vestigated  the  effects  of  preload  and  possible  damage  enhancement.  In  conjunc¬ 
tion  with  the  above  investigations,  fatigue  and  tensile  tests  are  being  per¬ 
formed  on  damaged  and  undamaged  specimens  to  establish  the  effect  of  impact  on 
these  properties. 

Damage  base  data  were  generated  using  nominal- thickness ,  fabricated  speci¬ 
mens  (the  leading-edge  thicknesses  are  identical  to  the  corresponding  blade 
type  at  the  50%  span  level).  Hal f- thickness  fabricated  specimens  were  then 
used  in  a  scaling-effects  investigation.  Preload  investigations  were  per¬ 
formed  on  the  fabricated  metal  specimens  to  simulate  the  effects  of  centri¬ 
fugal  loading  on  local  leading-edge  damage.  The  specimens  were  preloaded  at 
stresses  typical  of  the  impact  locations  (30  and  70%  span  levels). 


6.3  SUBTASK  C  -  MATERIAL-SCREENING  TESTS 
6.3.1  Introduction 


This  subtask  addressed  the  development  of  material-screening  tests  neces¬ 
sary  to  evaluate  and  effectively  rank  candidate  materials  for  fan  and  com¬ 
pressor  blades  that  possess  superior  F0D  resistance.  The  material  parameters 
investigated  in  this  study  include:  (1)  the  perforation  resistance  of  the 
material,  (2)  the  extent  the  material  may  plastically  deform,  (3)  the  extent 
to  which  the  material  is  vulnerable  to  catastrophic  structural  failure,  and 
(4)  the  extent  to  which  the  material  is  vulnerable  to  degradation  of  fatigue 
properties . 

In  this  investigation,  the  above  material  parameters  will  be  defined  in 
terms  of  response  to  impact  in  small-scale  tests.  The  perforation-resistance 
parameter  will  be  quantified  as  a  ballistic  limit  velocity,  VL,  for  a  given 
impact  configuration.  The  quantity  Vl  is  a  direct  measure  of  the  amount  of 
impact  energy  the  material  can  absorb  without  catastrophic  local  damage.  The 
plastic-deformation  parameter  will  be  quantified  by  conducting  impact  tests 
which  maximize  the  likelihood  of  plastic  deformation.  The  catastrophic  struc¬ 
tural-failure  parameter  will  be  quantified  in  a  Charpy  impact  test  designed  to 
produce  gross  structural  failure.  The  fatigue-properties  parameter  will  be 
quantified  in  terms  of  reduction  of  ultimate  fatigue  strength.  The  parameters 
defined  by  these  tests  may  be  considered  individually  or  combined  with  suit¬ 
able  weighting  functions  to  compute  an  overall  figure  of  merit  for  a  partic¬ 
ular  material . 
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6.3.2  Testing  and  Analysis 

The  experimental  work  in  this  subtask  has  been  completed,  and  the  data 
analysis  for  the  residual  fatigue  testing  is  underway. 

Table  10  gives  the  results  of  the  bal listic- limit  and  maximum- deforma¬ 
tion  testing.  The  projectiles  used  in  the  testing  included  the  0.25-inch 
diameter  granite  pebble,  the  0.5-inch  diameter  microballoon  gelatin  bird, 
and  the  0.5-inch  ice  ball. 


Table  10.  Ballistic  Limit  Results  (ft/sec). 


Target 

Mater ial 

0. 25-inch 
Pebble 

0. 5-inch 

Microbal loon 
Gelatin  Bird 

0 . 5-inch 
Ice  Ball 

410  Stainless 

3379 

>3012 

>2036 

Steel  (annealed) 

410  Stainless 

3335 

>1957 

— 

8-1-1  Titanium 

2859 

>1961 

>2173 

6-4  Titanium 

— 

>1923 

— 

B/Al  Composite 

633 

913 

G/Ep.  Composite 

— 

318 

The  specimen  size  used  in  the  testing  was  3x3  inch  having  a  thickness  of 
0.06  inches.  All  of  the  impacts  were  normal  impacts  (90°).  The  results  show 
that  the  composite  materials  had  the  lower  ballistic  limits. 

The  six  materials  evaluated  in  the  instrumented  Charpy  impact  tests  in¬ 
cluded  8-1-1  titanium,  6-4  titanium,  410  stainless  steel  in  the  annealed  con¬ 
dition,  410  stainless  steel  in  the  heat-treated  condition,  boron/ al uminum ,  and 
graphite/epoxy.  The  instrumented  Charpy  impact  tests  were  conducted  on  a 
Tinius  Olsen  pendulum  impact  machine  instrumented  with  a  Dynatype  loading  de¬ 
vice  and  associated  electronics.  Prior  to  testing,  dimensions  were  measured 
for  all  specimens.  Load-versus-time  and  energy-ver sus-time  curves  were  re¬ 
corded  on  a  storage  oscilloscope  and  photographed.  The  results  of  the  Charpy 
impact  test  analysis  are  shown  in  Table  11  for  individual  specimens.  The 
peak  force  (lb)  and  peak  energy  (ft- lb)  readings  were  measured  from  the  oscil¬ 
loscope  traces.  The  energy  must  be  corrected  to  account  for  the  reduction  in 
velocity  resulting  from  the  impact  with  the  sample.  The  last  column  of  Table 
11  gives  the  corrected  energy  (ft-lb)  taking  into  account  this  reduction  of 
velocity.  The  table  indicates  that  the  8-1-1  Ti  material  is  the  best  material 
in  regards  to  the  Charpy  test.  The  ranking  of  the  metals  would  be  as  follows: 
8-1-1  Ti  is  first,  followed  by  the  6-4  Ti  material,  then  the  heat-treated  410 
stainless  steel,  and  finally  the  410  stainless  steel  in  the  annealed  condi¬ 
tion.  For  the  composite  materials,  the  Charpy  tests  indicate  that  the 
graphite-epoxy  material  is  better  than  the  boron-aluminum  material. 
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Table  11.  Charpy  Test  Results. 


Hater ial 

Spec imen 
Number 

Pe  ak 

Force 

(lb) 

Peak 
Energy 
(ft- lb) 

Correc  ted 
Peak 
Energy 
(ft- lb) 

8-1-1  Ti 

1 

180.6 

10.1 

8.6 

8-1-1  Ti 

2 

183.6 

10.7 

9.1 

8-1-1  Ti 

3 

180.6 

10.7 

9. 1 

8-1-1  Ti 

4 

188.1 

10.9 

9.2 

6-4  Ti 

1 

16  7.2 

9.4 

8. 1 

6-4  Ti 

2 

16  7.2 

9.6 

8.3 

6-4  T  i 

3 

164.2 

9.0 

7.8 

6-4  Ti 

4 

165.7 

9.3 

8.1 

410  SS  (annealed) 

i 

77.6 

4.3 

4.0 

410  SS  (annealed) 

2 

79.1 

4.3 

4.0 

410  SS  (annealed) 

3 

76.1 

4.2 

3.9 

410  SS  (annealed) 

4 

76.1 

4.2 

3.9 

410  SS  (heat-treated) 

1 

153.7 

9.1 

7.9 

410  SS  (heat-treated) 

2 

153.7 

9.0 

7.8 

410  SS  (heat-treated) 

3 

149.3 

8.8 

7.7 

410  SS  (heat-treated) 

4 

149.3 

8.  7 

7.6 

Boron/ Aluminum 

C-  ID 

85 .1 

0.45 

0.45 

Boron/ Aluminum 

C-1E 

88.1 

0.51 

0.51 

Boron/ Aluminum 

C-2E 

92.5 

0.51 

0.51 

Boron/ Aluminum 

3D 

89.6 

0.51 

0.51 

Boron/ Alum  in urn 

C-3D 

87.3 

0.43 

0.43 

Graphite/Epoxy 

BA -3 -A 

148. 5 

0.87 

0.86 

Graphi te/Epoxy 

BA-3-B 

150.0 

0.97 

0.96 

Graphite/Epoxy 

BA-3 -D 

128.6 

0.  70 

0.69 

Graphite/Epoxy 

BA-3-G 

134.3 

0.78 

0.77 
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TASK  V  -  PARAMETRIC  ANALYSIS 


7.1  INTRODUCTION 


Parametric  relationships  describing  the  changes  in  dynamic  structural 
response  of  impacted  simple-plate  elements  with  the  progressive  introduction 
of  blade  airfoil  geometric  features  will  be  derived  in  this  task.  Addition¬ 
ally,  the  adequacy  of  the  structural-response  models  formulated  in  Tasks  II, 
III,  and  IV  will  be  investigated  and  modifications  made  as  necessary. 

This  task  is  broadly  divided  into  three  phases.  The  first  phase  is 
analysis  and  correlation  of  Task  VI  tests  conducted  with  simple,  cantilevered, 
flat  plates  with  different  aspect  ratios  and  thicknesses.  The  second  phase  is 
the  analysis  and  correlation  of  Task  VI  tests  with  geometry  and  impact  condi¬ 
tions  more  representative  of  actual  airfoil  FOD  dynamics  but  still  simple 
enough  that  the  sources  of  observed  effects  can  be  isolated.  The  third  and 
final  phase  of  this  task  is  the  analysis  and  correlation  of  Task  VI  and  Task 
IX  tests  of  the  three  selected  first-stage  airfoils  in  stationary  and  rotating 
environments . 


7.2  BLADES  SELECTED  FOR  ANALYSIS 


The  three  blades  selected  are  representative  of  first-stage  airfoils 
from  Air  Force  inventory,  production  development,  and  advanced-development 
engines.  These  blades  include  subsonic,  transonic,  and  supersonic  airfoil 
shapes  and  stainless  steel,  titanium,  and  an  advanced  composite  material.  The 
blades  selected  are  listed  below: 


Blade 

Represent ing 

Air f oi 1 

Material 

Shrouded 

379 

Air  Force  Inventory 

Subsonic 

Stainless  Steel 

No 

rioi 

Production  Development 

Transonic 

Titanium 

Yes 

APS  I 

Advanced  Development 

Supersonic 

B/Al  Composite 

No 

Supersonic  Engines 


7.3  PARAMETRIC  MATRIX 

A  parametric  matrix  has  been  formulated  to  define  the  impact  conditions 
and  response-model  geometries  that  will  be  analyzed  in  Tasks  V  and  X  and 
provide  guidance  in  the  selection  of  the  structural-element  and  full-scale- 
blade  tests  that  will  be  conducted  in  Tasks  VI,  IX,  and  XI. 
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7.4  NOSAPM  MODELING  AND  ANALYSIS 


7.4.1  Normal-Impact  Analysis 

The  following  paragraphs  discuss  the  modeling  and  analyses  of  Task  VI 
tests.  Table  12  lists  the  principal  conditions  of  one  of  these  tests.  It 
should  be  noted  that,  in  this  test,  both  the  projectile  mass  and  the  enclos¬ 
ing  sabot  impacted  the  blade.  Thus,  although  the  projectile  mass  alone  is 
3.81  cm  in  diameter  and  81.4  g  in  weight,  the  total  weight  of  100.5  g  with 
an  outside  diameter  of  5.08  cm  must  be  considered  for  full  momentum  transfer. 
The  geometry  of  the  test  plate  and  the  impacting  object  are  shown  in  Figure 
33.  The  planform  is  the  same  aspect  ratio  as  the  F101  blade. 

In  order  to  carry  out  the  analysis,  a  finite-element  model  of  this  test 
case  was  constructed  using  thirty-four  16-noded  elements.  The  finite-element 
model  is  shown  in  Figure  34  looking  straight  down  on  the  planform  of  the  plate. 
As  can  be  seen  in  the  figure,  the  mesh  was  made  finer  in  the  vicinity  of  the 
impact  in  order  to  better  approximate  the  expected  local  bulging.  Beaming 
action  is  expected  to  dominate  in  the  areas  remote  from  the  impact  site,  and 
the  elements  are  larger.  The  impact  pressure  is  spread  over  Elements  20,  21, 

24  and  25.  These  elements  bound  a  rectangular  area  of  16.94  cm2  (2.62  in.2). 
Because  the  sabot  as  well  as  the  projectile  mass  impacted  the  target  plate, 
there  is  a  degree  of  uncertainty  with  respect  to  the  details  of  the  loading 
condition.  It  was  decided  to  use  the  16.94  cm2  impact  area  which  is 
intermediate  to  the  total  area  of  mass  plus  sabot  and  mass  alone.  The  length 
of  the  impacting  mass  both  in  the  test  and  in  the  analysis  was  7.62  cm  (3 
in.).  The  properties  of  8-1-1  titanium  used  in  this  analysis  are  listed  on 
Figure  35  along  with  the  bilinear  approximation  to  the  stress-strain  curve 
used  for  the  finite-element  model.  It  should  be  pointed  out  that  these 
properties  are  more  appropriate  for  4-6  titanium  and  were  used  before  the 
8-1-1  data  from  UDRI  became  available.  As  a  result,  the  corrected  modulus  and 
yield  stress  should  be  used  in  future  analyses.  However,  these  properties  are 
adequate  for  gaining  modeling  experience. 

At  the  time  these  analyses  were  run,  the  interactive  loading  model  for 
normal  impact  was  not  yet  available.  As  a  result,  the  loading  history  applied 
to  Elements  20,  21,  24,  and  25  was  derived  from  simple  impulse-momentum  con¬ 
siderations  and  the  experimental  results  reported  by  Barber,  Wilbeck,  and 
Taylor6.  The  pressure  exerted  by  a  steady  jet  of  fluid  impacting  a  flat 
surface  can  be  calculated  as 

P  =  PV2  (21) 

where  p  is  the  mass  density  and  V  is  the  normal  impact  velocity.  Refer  to 
Figure  22. 

Since  the  test  under  consideration  involved  the  impact  of  both  sabot  and 
projectile  mass,  the  density  (p)  used  here  is  an  effect  '>e  density  calculated 
as  the  total  mass  of  sabot  and  projectile  (M)  divided  by  the  total  volume  of 
the  impacting  mass  (A  x  L).  Here  A  is  the  area  under  the  impact  on  the  model 
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(Elements  20,  2  L ,  24,  ami  25)  and  L  is  the  length  of  the  impactor.  Barber, 

W i  1  beck ,  and  Taylor  suggest  that,  during  a  normal  impact,  the  peak  pressure 
is  a  shock  pressure  about  twice  the  magnitude  of  the  steady-state  pressure  and 
occurs  at  roughly  20%  of  the  total  impact  time  duration.  Thereafter,  the 
pressure  reduces,  in  a  roughly  linear  manner,  to  zero;  refer  to  Figure  23. 
Using  this  approach,  the  model  for  the  impact  pressure  was  formulated.  This 
pressure  was  spread  uniformly  over  the  impact  area,  and  the  pressure  history 
is  exhibited  in  Figure  36. 

The  purpose  of  this  full-scale  analysis  was  to  develop  experience  in 
the  modeling  techniques  available  within  the  Level  1  package  for  application 
to  the  blade-impact  problem.  There  are  numerous  decisions  which  must  be  made 
in  conducting  such  a  finite-element  analysis,  including:  choice  of  an  appro¬ 
priate  time  step,  choice  of  the  number  of  Gauss  points  for  element-stiffness 
integration,  stiffness  reformation  strategy,  and  use  of  the  equilibrium 
iteration  technique  for  equation  solution.  A  great  deal  of  information  with 
respect  to  all  of  these  questions  was  gathered  during  the  course  of  the 
analyses . 


Table  12.  Dayton 

Shot  Number 
Date 

Projectile  Mass 

Total  Mass 

Projectile  Diameter 

Diameter  of  Projectile  and  Sabot 

Impact  Velocity 

Target  Material 

ImDact  Angle 

Center  of  Impact 
(as  percent  of  span  from  root) 

Mount 


test  Conditions. 

20094 

2-8 

81.4  g  (2.87  oz) 

100.5  g  (3.54  oz) 

3.81  cm  (1.5  in.) 

5.08  cm  (2.0  in.) 

177.4  m/sec  (582  ft /sec) 
8-1-1  Ti 

90° 

70% 

Cantilever;  31.12  cm  (12.25 
in . )  above  mount 
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figure  36.  Pressure-Loading  History 


Perhaps  Che  most  challenging  area  of  computation  strategy  in  this  pro¬ 
blem  is  the  extremely  wide  range  in  time  scales  that  must  be  considered 
during  a  complete  impact  analysis.  Table  13  lists  some  of  the  most  important 
characteristic  times  for  the  problem  described  herein  and  illustrates  the 
extreme  differences  in  time  frames.  As  can  be  seen,  these  time  frames  snan 
almost  four  orders  of  magnitude.  Yet  there  are  events  at  both  ends  of  the 
scale  which  are  regarded  as  significant.  For  example,  during  the  time  0  ♦  500 
ysec  the  impact  is  occurring,  and  the  intent  of  this  program  has  been  to  deal 
with  interaction  between  impactor  and  structure  during  this  time.  In  addition, 
it  is  during  this  time  that  large  local  stresses  appear  in  blades  and  plates 
and  cause  in  bulging  and  curlback.  On  the  other  hand,  the  large  tip  deflec¬ 
tions  and  large  root  stresses  associated  with  gross  structural  damage  occur  at 
the  other  end  of  the  time  scale:  10  to  20  msec  after  impact.  It  appears  that 
the  only  wav  to  dual  with  both  of  these  time  frames  is  with  two  time  steps. 

One  time  step  would  be  used  for  the  study  of  local  effects;  then,  using  the 
restart  capability  of  the  Level  l  package,  the  larger  time  frame  would  be 
studied  with  a  significantly  larger  time  step  using  the  last  step  of  the 
short-time  solution  as  initial  conditions. 


Table  13.  Characteristic  Times. 


Time  required  for  a  plane  stress  wave  to  travel  from 
centerline  to  plate  edge  (Y-face)  and  return  15  usee 


Time  required  for  a  plan  stress  wave  to  travel  from 
70"  span  to  tip  and  return  to  70%  span 

Characteristic  period  of  plate  response  for  the  square 
impact  area 

Time  required  for  a  plane  stress  wave  to  travel  from 
70%  span  to  root  and  return  to  70%  span 


40  usee 

80  usee 

90  usee 


Load  duration 

Characteristic  period  of  4th  flexural 
Characteristic  period  of  3rd  flexural 
Characteristic  period  of  2nd  flexural 
Time  to  maximum  tip  deflection 
Characteristic  period  of  1st  flexural 


430 

ps  ec 

mode 

870 

us  ec 

mode 

1740 

usee 

mode 

4900 

usee 

8000 

usee 

mode 

30,000 

usee 
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A  time  step  of  At  =  20  psec  was  chosen  for  this  analysis.  This  time  step 
provides  adequate  representation  for  the  loading  history  as  well  as  the  flex¬ 
ural  modes  of  the  problem.  It  is  too  large  to  closely  follow  the  movement  of 
plane  stress-wave  fronts  in  the  plate,  but  since  the  local  or  bulging  response 
is  nearly  quasi-static  with  respect  to  the  load  application  this  does  not 
appear  to  be  a  large  disadvantage.  For  this  analysis,  three  Gauss  points 
through  the  thickness  as  well  as  3  x  3  in  the  plane  of  the  plate  were  used  for 
numerical  integration.  There  is  some  indication  that  a  2  x  2  Gauss  point 
scheme  should  be  used  in-plane,  and  a  recommenation  on  this  issue  will  be  pur¬ 
sued.  If  effective,  it  would  substantially  reduce  the  cost  of  the  analysis. 
The  intervals  of  stiffness  reformation  and  the  use  of  equilibrium  iteration 
are  described  in  Figure  37.  Observation  of  the  change  in  geometry  as  well 
as  the  appearance  of  plasticity  in  the  problem  indicate  that  these  choices 
are  adequate.  It  has  been  discovered  during  application  of  NOSAPM  to  other 
structural  problems  at  CRD  that  the  equilibrium  iteration  technique  will  not 
work  during  large-scale  unloading  from  plastic  states  of  stress.  The  job  will 
prematurely  terminate  if  it  is  used  in  such  cases. 

Figures  38  through  42  illustrate  some  of  the  results  of  this  analysis 
for  times  less  than  500  psec  after  impact.  The  general  behavior  of  the 
response  as  well  as  the  nonlinear  nature  of  the  problem  are  illustrated. 

Figure  38  plots  the  lateral  displacement  of  the  plate  (X-displacement )  versus 
the  Y-location  at  the  70%  span  chord-line  designated  S-S  in  the  drawing.  The 
solid  lines  illustrate  the  full  nonlinear  response  of  the  plate  at  60,  100, 
220,  380,  and  500  psec.  As  can  be  seen,  at  60  and  100  psec  the  response  is 
characteristic  of  local  plate  bulging.  The  relatively  high-speed  flexural 
waves  associated  with  this  response  reach  the  edges  of  the  plate  around  80 
to  100  psec  after  the  impact;  thereafter  the  laterial  displacement  of  the 
plate  becomes  much  more  characteristic  of  a  beam  than  a  plate,  as  can  be  seen 
from  the  rather  constant  value  of  lateral  deflection  across  the  chord  line. 

The  dashed  lines  in  Figure  38  illustrate  the  linear  solution  to  the  same 
impact  problem.  At  60  and  100  psec  the  linear  and  nonlinear  solution  coin¬ 
cide.  However,  comparison  of  the  two  solutions  for  larger  times  indicates 
that  the  nonlinear,  in-plane  stiffness  associated  with  displacements  of  the 
order  of  the  plate  thickness  is  clearly  active  in  reducing  the  displace¬ 
ments  which  might  be  expected  from  a  linear  solution.  There  are  no  stresses 
above  the  yield  point  of  the  titanium  until  about  140  psec,  anu  even  then 
these  stresses  are  extremely  localized,  appearing  only  near  plate  surfaces 
where  tensile  membrane  and  bending  stresses  are  additive. 

Figure  39  plots  lateral  displacement  as  a  function  of  the  position 
along  the  plate  centerline  B-B  described  in  the  accompanying  figure.  The 
solid  lines  describe  the  nonlinear  plate  response  for  60,  140,  180,  300,  and 
500  psec  and,  again,  give  a  clear  indication  of  how  the  behavior  develops.  At 
60  psec,  the  response  is  characterized  by  local  bulging  as  previously  noted 
with  respect  to  Figure  38.  Between  about  140  and  180  psec  the  flexural 
response  reaches  the  blade  tip,  resulting  in  a  positve  tip  displacement.  By 
300  psec  some  of  the  high-frequency  flexural  waves  are  reaching  the  root  of 
the  plate.  With  reference  to  Figure  38,  it  should  be  noted  that  t.ie  plate 
response  is  now  much  more  characteristic  of  a  long  beam-plate  with  relatively 
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E  =  Equilibrium  Iteration 
R  =  Reform  Stiffness 
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Figure  37.  Stiffness  Reformation  and  Equilibrium  Iteration  Strategy. 
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Figure  38.  Lateral  Displacement  Alonn  Chord-Line. 
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constant  displacement  across  the  chord  line.  Finally,  by  500  jjsec  the  lower 
frequency  beaming  modes  are  beginning  to  develop,  and  the  tip  deflection  is 
now  negative.  For  comparison,  the  linear  solution  at  500  ysec  is  also  pre¬ 
sented  as  the  dashed  line.  As  pointed  out  with  reference  to  Figure  38,  the 
linear  solution  overestimates  the  flexural  displacements  of  the  plate  be¬ 
cause  of  the  absence  of  the  nonlinear,  in-plane  stiffness  developed  at  large 
displacements . 

Figures  40  and  41  help  emphasize  the  nonlinear  effects  present  in  the 
blade-impact  problem.  In  both  of  these  figures,  full  nonlinear  solutions  are 
compared  to  linear  solutions.  Figure  40  presents  the  lateral-deflection 
history  of  the  plate  at  the  very  center  of  impact.  Figure  41  presents  the 
tip-deflection  history.  In  both  cases  the  nonlinear  response  of  the  plate 
is  evident  in  the  form  of  additional  in-plane  stiffness  resulting  in  smaller 
displacements  with  respect  to  the  linear  solution.  The  effect  is  particularly 
evident  in  the  tip-displacement  history  of  Figure  41  were  the  nonlinear  effects 
result  in  a  maximum  positive  tip  displacement,  which  is  less  than  one-quarter 
that  predicted  by  linear  analysis.  As  indicated  in  Figure  41  plastic  stresses 
appear  after  about  140  psec,  but  they  remain  extremely  localized. 

Finally,  Figure  42  is  presented  to  illustrate  the  actual  shape  of  the 
deflected  plate  along  the  centerline  at  t  =  500  psec.  The  load  has  been  com¬ 
pletely  applied  by  this  time,  and  the  presence  of  the  first  t'.  -ee  cantilever- 
beam  modes  are  already  visibly  present  in  the  figure. 

The  analysis  was  refined  and  then  extended  from  the  time  frame  of  the 
load  application  to  include  the  longer  term  response  of  the  plate.  With 
respect  to  refining  the  previously  discussed  analysis,  two  items  should  be 
mentioned.  First,  a  10-psec  time  step  was  used  to  solve  the  problem  from  0  •* 
500  jsec  and  compared  to  the  solution  using  a  20-usec  time  step.  Sampling 
of  pertinent  results  indicated  that  the  reduction  in  the  time  step  decreased 
maximum  membrane  stresses  on  the  order  of  10%  while  increasing  the  maximum 
positive  tip  displacement  of  the  blade  (that  is  opposite  the  impact-loading 
direction)  by  approximately  30%.  The  spatial  distribution  of  membrane  stresses 
under  the  load  was  also  smoothed  out  using  the  smaller  time  step.  In  addition 
to  this  comparison,  the  number  of  Gauss  points  used  for  numerical  integration 
in  the  plane  of  the  plate  was  also  altered  from  a  3  x  3  configuration  to  a 
2x2  conf igviration  for  the  sake  of  comparison.  While  the  3x3  set  of  Gauss 
points  resulted  in  extreme  discontinuities  in  stresses  between  elements,  the 
2x2  scheme  provided  very  smooth  and  well-behaved  distributions.  Zienkiewicz 
also  notes  in  his  text^  that  lower  order  forms  of  Gaussian  integration 
may  sometimes  be  more  appropriate  in  terms  of  stress  distributions.  As  a 
result,  it  is  recommended  that  a  2  x  2  scheme  of  Gaussian  integration  be  used 
in  the  plane  of  the  plate  and  three  Gauss  points  through  the  thickness. 

An  indication  of  some  of  the  stress-response  behavior  observed  in  this 
analysis  is  displayed  in  Figures  43  and  44.  These  plots  focus  upon  the  non¬ 
linear  membrane  stresses  generated  in  the  plate  during  impact.  These  are 
effects  that  would  not  be  predictable  with  a  linear  analysis.  Figure  43 
illustrates  the  azz  membrane  stress  generated  near  the  center  of  impact 
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Figure  43.  Membrane  Stress  History  at  Center  of  Impact 


Membrane  Stress 


:i  a  function  of  time.  As  can  be  seen,  before  approximately  60  usee  there 
is  no  membrane  stress  response  since  the  displacements  in  the  area  of  the 
impact  are  less  than  the  plate  thickness,  and  linear  bending  theory  applies. 
However,  shortly  thereafter,  as  the  displacement  continue  to  grow,  a  large 
tensile  membrane  stress  begins  to  develop  as  a  result  of  nonlinear  effects. 
This  stress  passes  through  a  maximum  of  over  50  ksi  before  decaying  and  be¬ 
ginning  to  reflect  the  higher  frequency  response  associated  with  reflections 
from  the  plate  end  and  root.  This  pulse  in  the  membrane  stress  spreads 
and  attenuates  as  it  propogates  down  the  length  of  the  plate.  By  the  time  it 
reaches  the  root  area,  the  magnitude  is  significantly  reduced;  this  can  be 
seen  in  Figure  44  which  plots  the  azz  membrane  stress  history  at  the  Gauss 
point  nearest  to  the  centerline  at  the  root.  The  maximum  stress  here  is 
approximately  13  ksi.  Figure  45  plots  the  centerline  displacements  of  the 
plate  at  three  different  time  intervals  where  plastic  action  was  occurring. 

At  0.320  msec  the  yielding  is  local  to  the  impact  area.  At  1.15  msec  yield¬ 
ing  begins  to  occur  at  the  root.  This  yielding  seems  to  be  associated  with 
the  arrival  of  a  flexural  stress  wave  associated  with  the  third  natural 
bending  frequency  of  the  plate.  Unloading  then  occurs  at  the  root,  and  yield¬ 
ing  eventually  reappears  this  area  near  the  maximum  tip  displacement 
reached  at  approximately  6.2  msec. 


7.4.2  Comparison  Between  First-  and  Second-Level 
Analyses  for  Normal  Impact 


Further  refinement  relative  to  the  material  parameters  was  accomplished, 
and  a  comparison  between  NOSAPM  and  COMET  (second-level  analysis)  was  made  for 
normal  impact.  Material  data,  based  upon  tensile  tests  carried  out  at  UDRI  at 
various  strain  rates  in  Task  IV-A,  were  used  in  the  anlayses.  previous  ex¬ 
perience  indicates  that  strain  rates  in  these  impacted  plates  will  be  in  the 
range  of  375  in. /in. /sec.  With  this  prediction  and  the  UDRI  test  data,  the 
following  material  properties  were  used  in  the  analysis: 

=  17.3  x  10^  nsi 

E2  =  0.376  x  10^  psi 

v  =  0.3 

ay  =  191  ksi 

p  =  0.158  lb/in.3 


In  this  case  is  the  linear  elastic  Young's  modulus  while  Et  is  the 
tangent  modulus  which  approximates  the  material  response  after  yielding  has 
occurred.  The  same  NOSAPM  finite-element  model  was  used  for  this  analysis  as 
in  the  previously  described  normal-impact  analyses.  The  pressure  over  the 
impact  area  was  calculated  to  be  7100  psi  and  occured  90  usee  after  the 
beginning  of  imr’.ct.  The  total  impact  event  was  assumed  to  last  430  usec  . 
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For  the  Level  1  analysis,  a  10-ysec  time  step  was  used.  For  the  Level 
2  analysis,  a  l-psec  time  step  was  used.  The  two  results  agree  quite  closely 
over  the  time  frame  examined  in  spite  of  the  fact  that  the  elements  and  nu¬ 
merical  techniques  are  quite  different.  Unfortunately,  the  Level  2  analysis 
did  not  prove  to  be  substantially  less  expensive  than  the  Level  1  analysis  for 
the  computing  strategies  incorporated.  In  an  effort  to  further  test  the 
choice  of  time  step  used  in  the  Level  1  anlaysis,  the  same  problem  was  studied 
over  a  reduced  time  frame  wi.h  a  1-psec  time  step.  At  this  time  step  there 
was  virtually  no  visible  difference  between  the  10-psec  and  1-psec  solutions. 
Perhaps  a  more  pertinent  convergence  test  in  this  case  is  a  comparison  of 
the  membrane  stress  m  the  impact  area  since  the  propagation  of  this  stress  is 
at  near  the  speed  of  sound  in  the  plate  and,  hence,  takes  place  very  quickly. 
Over  the  time  frame  studied,  the  maximum  membrane  stress  under  the  load  was  on 
the  order  of  6L  different  for  the  1-psec  time-step  strategy  compared  to  the 
10-psec  investigation.  Thus,  although  improvements  in  the  membrane-stress 
distributions  can  be  made  by  reducing  the  time  from  10  psec  to  1  ysec,  one 
must  trade  off  this  benefit  against  cost  considerations. 


7.4.3  Comparison  Between  First-Level  Analysis  and  Experimental 
Data  for  Normal  Impact 


Two  normal-impact  correlations  were  carried  out  with  experimental  data 
obtained  by  UDRI  in  Tasks  III  and  VI.  In  both  cases  the  impact  angle  was 
90°;  as  a  result,  a  noninteractive  impact  pressure  characterization  (based 
upon  previously  reported  principles)  was  applied.  The  correlation  obtained 
in  these  two  examples  indicates  the  success  of  NOSAPM  in  predicting  dis¬ 
placement  histories  and  deformed  shapes  both  during  the  initial  phases  of 
impact  when  local  deformation  governs  as  well  as  time  periods  extending  to 
maximum  tip-displacement  intervals.  It  must  be  noted,  however,  that  achiev¬ 
ing  maximum-tip-displacement  results  is  much  more  expensive. 

The  geometries  of  the  two  normal-impact  problems  under  consideration  here 
are  characterized  in  Figure  46.  All  dimensions  associated  with  the  two  tests 
are  identical  except  for  the  plate  thickness;  that  is  4.27  mm  in  the  first 
case  and  2.21  mm  in  the  second  case.  The  impact  masses  and  velocities  are 
also  different.  The  first  impact  discussed  involved  an  impact  mass  of  100  g 
and  a  velocity  of  177.4  m/sec.  The  second  test  used  a  mass  of  82.6  g  impact¬ 
ing  at  103.8  m/sec.  The  finite-element  model  illustrated  in  Figure  47  has 
been  discussed  previously.  The  impact  area  in  this  model  is  comprised  of 
Elements  20,  21,  24,  and  25.  Based  upon  University  of  Dayton  Research  Insti¬ 
tute  data,  the  Young's  -.odulus  for  the  8-1-1  titanium  used  in  this  experiment 
was  shown  to  be  119  GPa  (17.3  x  10^  psi)  with  a  secondary  modulus  of  2.59  GPa 
(385  ksi).  The  material  density  used  in  the  analysis  was  4.38  g/cm-'  (0.158 
lb/in.  ^).  Although  strain  rates  of  200  to  300  sec  ^  have  been  measured 
and  predicted  analytically  during  impact  response,  preliminary  analyses 
indicated  that  the  most  substantial  yielding  in  these  tests  occurs  at  strain 
rates  of  the  order  of  4  sec-'-.  As  a  result,  the  1.13  GPa  (164  ksi)  yield 
stress  applied  in  this  analysis  reflects  such  a  strain  rate.  The  results 
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presented  in  Figures  48  and  49  were  obtained  with  a  2  x  2  Gauss  point  array 
in  the  plane  of  the  plate  and  four  Gauss  points  through  the  thickness. 

Stagliano  and  Mente^  feel  that  four  Gauss  points  are  necessary  through  a 
thin  plate  in  order  to  accurately  describe  plastic  deformation  incurred  during 
large-displacement  bending.  Results  from  our  work  using  both  three  and  four 
Gauss  points  through  the  thickness  are  inconclusive  in  this  respect. 

Figure  48  compares  the  history  of  the  tip-deflection  prediction  analytic¬ 
ally  with  the  data  collected  using  high-speed  photography  during  the  actual 
test.  As  can  be  seen,  the  correlation  is  reasonably  good,  especially  when  one 
considers  the  extreme  nonlinearity  of  the  problem.  The  analysis  predicts 
larger  tip  deflections  than  those  actually  measured  during  the  time  frame 
between  1  and  3  msec  after  impact  initiation.  It  underpredicts  the  maximum 
deflection  by  only  about  13%.  The  discrepancy  is  most  likely  a  reflection 
of  the  approximations  made  in  describing  the  load  history  and  the  "discretiza¬ 
tion"  applied  to  model  the  plate. 

During  the  second  normal-impact  test  considered  here  (Shot  40055),  at¬ 
tention  was  focused  upon  the  local  response  of  the  plate.  Moire  fringe 
techniques  were  used  to  get  a  detailed  description  of  the  local  deformation 
during  the  impact  response  (refer  to  Figure  25).  The  plate  impacted  in  this 
test  had  the  same  planform  dimensions  but  was  only  2.21  mm  thick.  The  pro¬ 
jectile  mass  was  82.6  g,  and  the  impact  velocity  was  103.8  m/sec.  Application 
of  the  same  principles  previously  discussed  with  respect  to  modeling  the  im¬ 
pact  loading  results  in  a  peak  pressure  of  14.0  MPa  at  150  usee  and  a  total 
impact  load  duration  of  730  jjsec.  Figure  49  compares  the  results  of  this  test 
with  analytical  predictions.  The  displacement  of  the  centerline  of  the  plate 
is  plotted  for  two  different  times  after  impact  initation:  150  and  460  usee. 
The  analysis  overestimates  the  local  maximum  of  the  plate  response  by  about  30% 
at  460  ysec  after  the  impact  initiation.  The  most  likely  reason  for  this  over¬ 
estimation  is  inaccuracy  in  the  load  model.  The  simple  approach  used  herein 
to  obtain  impact  pressures  is  based  upon  the  assumption  that  the  target  is 
rigid.  In  the  case  of  a  thin  plate  this  is,  of  course,  not  true.  The  plate 
develops  substantial  structural  velocities  during  an  impact,  and  this  response 
will  reduce  peak  pressure  loads  and  increase  the  duration  of  the  impact.  In 
addition,  although  the  impact  areas  over  which  the  pressure  loads  were  applied 
in  the  analysis  was  roughly  the  cross-sectional  area  of  the  projectile,  it  has 
been  observed^*  that,  in  a  normal  impact,  there  are  substantial  pressures  over 
an  area  as  much  as  four  times  the  cross  section  of  the  projectile.  Proper 
accounting  of  both  of  these  effects  would  tend  to  reduce  the  predicted  maximum 
local  displacement. 


7.4.4  Oblique-Impact  Analysis 

The  following  discussion  pertains  to  the  results  obtained  from  the  com¬ 
bined  program  representing  the  UDRI  loading  model  and  NOSAPM  response  model. 
An  example  of  the  pressure  calculation  for  a  24.4*  impact  on  a  flat  plate 
with  a  bird  velocity  of  15,240  in. /sec  is  presented  in  Figure  50.  The 
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Figure  48.  Tip-Displacement  History:  NOSAPM  Predictions  Vs. 
Test  for  Simulated-Bird  Impact. 
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Figure  50.  Pressure  Distribution  for  9  =  24.4°  Impact  on  a  Flat  Plate. 
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density  of  the  substitute-bird  material  used  in  the  analysis  was  0.03567  lbm/ 
in. ^  The  computed  pressure  distributions  were  symmetric  about  the  centerline 
in  the  span  direction  during  the  initiation  of  impact,  as  would  be  expected. 

The  maximum  response  occurred  at  the  leading  edge  and  decayed  to  zero  to¬ 
ward  the  trailing  edge.  With  35  pressure  nodes,  the  maximum  pressure  at  the 
leading  edge  was  7217  psi,  almost  exactly  twice  the  momentum-balance  pressure: 


(p/g)V2  sin2  8  =  [0. 03567/32. 2( 12) ] ( 15240)2(sin  24. 4*)2  =  3659  psi 

.22) 


This,  in  turn,  means  that  the  momentum  transfer  is  correct,  at  least  at  time 
zero.  The  predicted  pressures  are  somewhat  lower  if  a  coarser  grid  is  used 
to  define  the  pressure  distribution.  Figure  51  shows  the  impact  geometry 
and  pressure  distribution  and  Figure  50  shows  the  details  of  the  pressure 
distributions  for  two  pressure-grid  sizes. 

Solution  checks  were  also  run  on  the  following  operations;  in  all  cases, 
the  results  were  positive. 

•  Operation  of  load  interface  with  centrifugal  loads. 

•  Validity  of  calculation  of  equivalent  nodal  loads  from  updated 
pressure  data. 

•  Updating  of  pressure  data  during  impact. 

•  Termination  of  impact  loading. 

7.5  COMET  MODELING  AND  OBLIQUE-IMPACT  ANALYSIS 


The  loading  model  developed  by  CR&DC  (described  in  Section  5.4.5  of  this 
report)  has  been  incorporated  into  the  COMET  response  program,  and  initial 
checks  have  been  made.  This  loading  model  takes  account  of  shock  waves  and 
other  mechanisms,  described  in  Reference  6,  while  treating  the  bird/blade 
relationship  interactively.  The  approach  automatically  takes  account  of 
changes  in  blade  angle  during  impact  and  of  lengthening  of  the  impact 
duration  by  the  blade  displacement. 

Results  using  this  routine  are  presented  in  Figures  52  through  55  for 
a  3/8-in.  thick  aluminum  blade  model  extending  from  R  *  10  in.  to  R  =  20  in. 
The  blade  chord  extends  5  in.  The  bird  length  is  3  in.,  and  the  impact  region 
extends  from  R  *  15.5  in.  to  R  *  18.5  in.  The  analysis  was  run  using  66  ele¬ 
ments  and  66  nodes  at  1-in.  spacing  with  a  time  step  of  0.000002  sec  for  200 
times  steps. 
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Figure  52.  Model  of  Blade/Bird 


Impact  for  COMET  Analysi 


Lateral  deflections  at  the  leading  edge  at  0.0004  seconds  reached  about 
an  inch  in  the  axial  and  the  circumferential  directions.  The  lateral  deflec¬ 
tions  were  maximum  in  the  impact  region  and  were  lower  at  the  tip. 

The  deflection  at  the  17-in.  chord  line  show  that  initially  the  centri¬ 
fugal  forces  rotate  the  blade  a  small  amount  into  the  plane  of  rotation.  Sub¬ 
sequently  the  bird-impact  forces  bend  the  blade  back  substantially. 

The  interaction  of  bird  and  blade  is  shown  in  Figure  55.  It  is  evident 
that  blade  deflection  extends  the  duration  of  impact  substantially.  The  nomi¬ 
nal  impact  time  is  0.000278  sec  (3  in.  bird  length  and  10,786  in. /sec  bird 
relative  velocity),  but  the  figure  shows  the  impact  is  still  in  progress  at 
0.0004  sec  due  to  blade  flexing. 

The  peculiar-shape  bird  (highly  angled)  shown  in  the  figure  is  the  result 
of  a  need  by  the  University  of  Dayton  to  model  the  front  of  the  bird  parallel 
to  the  blade.  This  has  been  followed  in  the  present  COMET  version;  however, 
the  programming  in  COMET  can  permit  the  front  and  back  of  the  bird  to  have 
other  configurations. 
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8.0  TASK  VI  -  STRUCTURAL  ELEMENT  TESTS 


8.  1  INTRODUCTION 

The  work  of  Task  VI,  "Structural  Element  Tests,"  is  concerned  with  per¬ 
forming  nonrotating  bench  impact  tests  on  specimens  ranging  from  simple 
elements  to  real  blades.  The  simple  elements,  such  as  beams  and  plates, 
will  be  tested  with  progressive  introduction  of  airfoil  geometric  parameters 
to  validate  experimentally  the  analytical  predictions  of  Tasks  V  and  VIII. 
The  purpose  of  Task  V  is  to  derive  parametric  relationships  describing  the 
changes  in  dynamic  structural  response  as  a  function  of  material  and  geo¬ 
metric  features.  The  purpose  of  Task  VIII  is  to  derive  criteria  for  predic¬ 
ting  blade  FOD  tolerance.  These  criteria  will  be  formulated  to  make  full 
use  of  the  transient  dynamic  analysis  and  experimental  results  obtained  from 
other  tasks  of  the  program. 


8.2  DESCRIPTION  OF  TESTS 

The  materials  being  investigated  in  the  structural  element  tests  include 
403  stainless  steel,  8-1-1  titanium,  and  the  boron/ aluminum  composite.  The 
geometry  of  the  test  specimens  is  arbitrary  for  the  purpose  of  checking  out 
the  analytical  predictions  of  Tasks  V  and  VIII;  however,  the  geometry  and 
thickness  were  selected  to  be  typical  of  actual  blades  of  the  same  material, 
at  the  50%  span  location  along  the  leading  edge. 

The  impactors  for  the  study  include  2-  to  3-ounce  birds,  1.5-pound 
birds,  and  3-inch  ice  cylinders.  Microballoon  gelatin  material  is  used  to 
simulate  real  birds.  The  impact  velocities  were  selected  to  cause  no  damage, 
threshold  damage,  and  severe  damage.  The  impact  angles  selected  were 
arbitrary;  however,  they  will  correspond  to  those  of  actual  blades  at  the 
30  and  70%  span  locations. 

The  blade  parameters  being  investigated  include  aspect  ratio,  thickness/ 
chord  ratio,  shape  (constant  thickness,  wedge,  and  blade-like),  camber,  twist, 
and  shrouds. 


8.2.1  Composite  Structural  Specimens 

The  composite  specimens  include  flat  panels,  airfoil-like  cross  sections, 
cambered  shapes,  cambered  and  twisted  shapes,  and  real  APSI  blades.  The  fab¬ 
ricated  composite  specimens  are  shown  in  Figure  56. 


8.2.2  Metal  Structural  Specimens 

The  metal  specimens  include  flat  plates,  airfoil-like  cross  sections, 
cambered  shapes,  cambered  and  twisted  shapes,  and  real  J79  and  F101  blades. 


V. 
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Cross-Ply  B/Al 
(0/22/0/-22) 
Configuration  (1) 

Configuration  (2) 

Configuration  (3) 


Blade  Type  Aspect  Ratio 
Panel 

Same  as  (1) ,  Except  Length 
Reduced  by  1/2 
Same  as  (1),  Except 
Thickness  Reduced  by  1/2 


(a)  Configurations  1,  2,  and  3 


0.155  in. 


7.532 

in. 


J79  Specimen  Die  Used  and 
Aspect  Ratio  Same  as  APSI 
Configuration  (6)  -  Cambered  Panel 
Configuration  (7)  -  Same  as  (6), 
Except  Twist 
Introduced  After 
Forming 


(c)  Configurations  6  and  7 


Figure  56.  Task  VI  B/Al  Composite  Test  Specimen  Configurations. 
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The  titanium  material  was  selected  for  a  baseline  series  of  tests,  and 
the  stainless  steel  was  selected  for  a  supplementary  series  of  tests.  The 
baseline  series  of  tests  on  the  fabricated  titanium  specimens  considered  all 
geometrical  effects  except  camber  and  twist.  It  has  been  established  that 
camber  and  twist  are  not  cold-formable  on  the  titanium  specimens  and  that 
these  parameters  could  be  incorporated  only  through  hot-forming  or  machining 
operations.  Therefore,  camber  and  twist  were  investigated  utilizing  the 
stainless  steel  cold-formed  specimens. 


8. 3  DISCUSSION 

All  of  the  impact  tests  for  Task  VI  have  been  completed,  and  analyses 
of  the  strain  data  and  of  the  displacement  responses  recorded  on  high-speed 
movies  are  underway.  Table  14  describes  some  of  the  tests  performed  on  real 
blades.  During  the  testing  of  the  fabricated  structural  specimens,  the 
different  geometry  effects  (i.e.  aspect  ratio,  thickness/chord  ratio,  etc.) 
were  introduced  in  a  different  sequence  for  the  titanium  and  steel  specimens 
so  that  any  synergistic  effects  that  might  exist  could  be  accurately  assessed. 

The  strain  data  will  be  formated  to  obtain  strain-  and  strain-rate-ver sus- 
time  plots  for  direct  comparison  with  the  predictions. 

Detailed  analysis  of  the  high-speed  movies  is  being  conducted  to  obtain 
the  d isplacement-versus-time  profiles  for  comparison  with  the  predicted  dis¬ 
placement  response. 

Figures  57  through  59  show  examples  of  the  displacement,  strain,  and 
strain-rate  responses  that  were  obtained  from  a  3-ounce-bird,  normal- impac t 
test  at  the  70%  span  location  on  a  8-1-1  Ti  plate-shaped  specimen  with  a 
blade-like  (F101)  aspect  ratio. 
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Table  14.  Some  Task  VI  Impact  Tests  Conducted  on  Real  Blades. 


'0.00  20.00  40.00  60.00  BO. 00  100.00 

Time,  msec 


Figure  58.  Strain/Time  Curve  for  Root  Gage 
in  Span  Direction 


9.0  TASK  VII  -  ERROR-BAND  ANALYSIS 


9.1  INTRODUCTION 


The  objective  of  this  task  is  to  analytically  determine  the  uncertainties 
in  strain  history,  deflection  history,  and  damage-related  parameters  of  FOD 
which  are  introduced  by  such  factors  as  manufacturing  tolerance,  variability 
in  material  properties,  and  variability  in  impact  conditions.  This  determina¬ 
tion  is  important  in  Task  VIII  for  quantifying  variances  for  the  damage  model 
and  in  Task  VI  and  IX  for  comparing  computer  output  with  test  results. 

The  damage-related  parameters  will  be  among  those  used  in  Task  VIII  as 
the  basis  for  the  impact  damage  model  (maximum  stress  or  maximum  strain  for 
example).  These  parameters  will  be  obtained  as  output  from  the  analysis  models 
defined  in  Task  II.  It  is  expected  that  the  second-level  response  model  will 
be  adequate,  but  the  detailed  modeling  required  for  studying  the  effects  of 
local  deviations  in  geometry  within  manufacturing  tolerance  might  require  the 
use  of  the  first-level-analysis  reponse  model. 


9.2  PROPERTIES  REQUIRING  AN  ERROR-BAND  ESTIMATION 

A  list  of  geometries  has  been  compiled  relating  to  blade  material,  impac- 
tor  characteristics,  and  aerogeometric  factors  for  which  the  FOD  design  cri¬ 
teria  (developed  in  Task  VIII)  will  most  need  an  error-band  estimation.  It  is 
believed  that  variation  of  these  properties  within  a  reasonable  range  can 
cause  a  wide  variation  in  FOD.  These  error-band  properties  and  the  associated 
damage  modes  are  identified  in  the  matrix  shown  on  Table  15. 


9 . 3  DAMAGE  VECTOR 


The  concept  of  a  damage  vector  will  be  utilized  in  the  response  models  to 
quantify  the  damage  modes  identified  previously.  This  damage  vector  will  be 
used  in  conjunction  with  the  sensitivity  studies  for  the  error-band  properties. 
The  damage  vector  will  consist  of  the  numerical  values  of  perhaps  a  dozen  quan¬ 
tities  which  signal  the  onset  of  each  of  the  damage  modes.  For  example,  rela¬ 
tive  to  blade  displacement  and  local  tearing,  the  locations  and  values  of  the 
maximum  and  minimum  axial  and  tangential  blade  displacements  and  the  location 
and  values  of  the  maximum  principal  strains  will  be  included  in  the  damage 
vector . 

The  damage  vector  is  essentially  a  list  of  the  values  of  the  failure  cri¬ 
teria.  This  vector  provides  the  user  with  an  appraisal  of  the  damage  sus¬ 
tained.  The  degree  to  which  the  damage  vector  is  altered  when  changes  are 
made  in  material,  impactor,  or  aerogeometric  parameters  would  identify,  in  an 
unambiguous  way,  the  sensitivity  of  each  damage  mode  to  changes  in  these 
parameters . 


Table  15.  Matrix  of  Error-Band  Properties  and  Associated  Damage  Modes. 
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LO.O  TASK  VIII  -  FOREIGN  OBJECT  IMPACT  DESIGN  CRITERIA 


10. I  INTRODUCTION 

-  I 

The  purpose  of  this  task  is  to  derive  the  criteria  for  predicting  blade 
FOD  tolerance.  The  criteria  will  be  formulated  to  make  full  use  of  the  tran¬ 
sient  dynamic  analysis  and  experimental  results  obtained  in  the  previous  tasks 
and  will  be  based  on  a  combination  of  theory  and  regression  analysis.  Dimen¬ 
sional  analysis  and  an  examination  of  the  physics  of  impact  will  determine  i 

which  impact-dynamics  and  airfoil-response  parameters  will  be  used  to  formu¬ 
late  the  damage  model.  Regression  analysis  will  be  used  to  determine  the 
constants  in  the  model  by  correlation  with  the  experimental  results.  The 
experience  gained  at  General  Electric  in  the  development  of  simplified  cri¬ 
teria  based  on  purely  elastic  analyses  will  provide  an  excellent  basis  for 
the  complete  elastic-plastic  analysis  of  this  program. 

The  damage  models  formulated  will  be  correlated  with  the  observed  damage 
from  Tasks  VI  and  IX  tests,  and  each  will  be  sufficiently  resolved  to  predict 
a  threshold  value  beneath  which  the  damage  is  not  observed  and  above  which 
the  prescribed  damage  is  evident.  The  level  of  damage  for  such  definitions 
as  local  plastic  deformation  will  be  of  relatively  increasing  severity  from 
no  damage  to  large  plastic  deformation  so  that  multiple  qualitative  thresholds 
will  be  required.  A  more  definitive  damage  criterion  is  to  specify  the  size 
of  the  region  of  local  plastic  deformation  relative  to  a  normalizing  dimension 
of  the  airfoil  -  such  as  the  chord.  Damage  definitions,  such  as  local  frac¬ 
ture,  are  qualitatively  more  precise  measurements;  therefore,  a  single  thresh¬ 
old  will  be  sufficient. 


10.2  DESIGN  CRITERIA 


An  assessment  has  been  completed  of  how  well  the  structural-response 
models  can  provide  the  selected  failure  variables.  It  was  established  that 
the  NOSAPM  program  in  the  current  form  provides  stresses  only  at  the  Gauss 
points  for  nonlinear  problems.  Because  a  forming-limit  concept  has  been 
selected  that  uses  the  relative  values  of  the  maximum  and  minimum  principal 
strains  to  establish  whether  tearing  occurs  in  metal  blades,  and  also  because 
the  failure  criteria  selected  for  composite  blades  is  based  on  the  maximum 
bending  stresses  and  shear  stresses,  the  NOSAPM  program  will  be  modified  to 
provide  surface  principal  strains  and  physical  stresses  as  well  as  the  Gauss- 
point  stresses.  In  addition,  surface  physical  strains  will  be  provided  to 
allow  correlation  with  strain  gage  data  obtained  from  the  Task  VI  impact 
tests.  The  COMET  program  provides  physical  and  principal  surface  strains  and 
surface  stresses. 

Both  NOSAPM  and  COMET  provide  geometric  and  material  nonlinear  analysis 
capability  that  addresses  the  large-displacement,  e lastopl as t ic  behavior  that 
is  characteristic  of  the  transient  response  of  first-stage  fan  and  compressor 
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blades  due  Co  foreign-object  impact.  This  geometric  and  material  nonlinear- 
analysis  capability  will  be  used  to  predict  local  plastic  deformation,  such 
as  bulging,  and  overall  blade  displacements  resulting  from  far-field  struc¬ 
tural  response. 

The  design  criteria  that  have  been  selected  will  be  consistent  with  the 
predictive  capabilities  of  the  Task  II  transient  structural-response  models 
and  with  the  test  data  and  methods  of  Task  IV.  They  are  presented  below. 

10.2.1  Design  Criteria  for  Use  with  Task  II  Transient  Structural-Response 
Models  and  Task  IV-A  Test  Data 

Metal  Blades 

«  Limit  local  plastic  deformation  to  a  specific  percentage  of  chord 
to  avoid  stall/surge. 

•  Limit  strains  in  the  impact  region  in  order  to  avoid  tearing  and 
mass  loss. 

•  Limit  strains  at  the  root  of  the  blade  in  order  to  prevent  airfoil 
loss . 

•  Limit  displacements  in  order  to  avoid  contact  with  adjacent  blades 
and  stages. 

Composite  Blades 

•  Limit  spanwise  and  chordwise  normal  stress  in  order  to  avoid  fiber 
failure. 

•  Limit  shear  and  through-the-th icknes s  normal  stresses  in  ord^r  to 
avoid  delamination. 

•  Limit  displacements  in  order  to  avoid  contact  with  adjacent  blades 
and  stages. 

•  Limit  spanwise  normal  stresses  at  the  root  of  the  blade  in  order 
to  prevent  airfoil  loss. 

10.2.2  Design  Criteria  Based  on  Task  IV-A,  B,  and  C  Test  Data  and  Methods 

•  Utilize  bank  of  design  data  to  relate  residual  tensile  and  fatigue 
properties  to  local,  leading-edge,  ballistic-impact  damage. 

•  Evaluate  and  rank  candidate  blade  materials  through  screening  tests 
that  address  the  ballistic  limit,  local  deformation  charact eris t ics , 
gross  structuraL-damage  characteristics ,  and  residual  fatigue 
strength. 
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11.0  TASK  IX  -  SINGLE-BLADE  IMPACT  TESTS 


11.1  INTRODUCTION 


The  effort  of  Task  IX,  "Single  Blade  Impact  Tests,"  is  concerned  with 
corroborating  analytical  and  parametric  predictions  and  evaluating  the  valid¬ 
ity  of  design  criteria  developed  in  other  tasks.  This  will  be  accomplished 
by  impact  testing  full-scale  instrumented  blades  in  an  existing  whirligig 
facility. 

Three  blade  designs  selected  for  transient-response  impact  analysis  in 
Task  V  will  be  tested  under  test  conditions  that  duplicate  the  analytical 
impact-parameter  values  used  in  the  earlier  tasks.  This  is  necessary  to 
obtain  a  direct  correlation  of  the  analytical  and  experimental  results. 

Activities  conducted  under  this  task  will  be  carefully  coordinated  witn 
the  analytical  studies  of  Task  II,  Task  V,  and  Task  VIII  and  the  experimental 
evaluation  of  Task  VI. 


11.2  TEST  PLAN 


impact  testing  of  full-scale,  instrumented  blades  will  be  conducted  in 
the  existing  whirligig  facilities.  The  three  blade  designs  selected  in  Task  V 
will  be  tested.  These  are  the  APSI  fan  B/Al  blade,  the  J79  compressor  Stage  1 
steel  blade,  and  the  F101  fan  Stage  1  tip-shrouded  titanium  blade.  The  J79* 
and  APSI  blades  will  each  be  tested  using  an  APSI/J79  single-blade  test  disk 
set  up  in  General  Electric,  Evendale,  test  cell  No.  63.  The  F101  blades  will 
be  tested  using  an  FlOi  engine  disk  setup  with  a  soft-mount  TF34  vehicle  in 
General  Electric,  Evendale,  test  cell  No.  62.  Figure  60  shows  the  single¬ 
blade  impact  setup  for  the  APSI  and  J79  blades,  and  Figure  61  shows  the  setup 
for  the  F101  blade. 

Four  blades  of  each  of  the  three  selected  configurations  will  be  impact 
tested.  Each  blade  will  be  set  up  for  a  single-blade  impact  test  and  instru¬ 
mented  with  three  dynamic  strain  gages. 

Rotational  speed  and  angle  for  the  J79  blades  will  be  set  to  simulate  a 
bird  ingestion  at  takeoff  conditions.  A  simulated  bird  (or  an  ice  slab)  will 
be  injected  into  the  path  of  the  rotating  blade.  An  existing,  angled  slot  in 
the  disk  provides  additional  blade-twist  angle  to  compensate  for  the  zero- 
velocity  state  of  the  projectile  on  impact  and  to  provide  the  proper  blade-to- 
projectile  incidence  angle. 

Since  it  is  planned  to  use  existing  test  hardware  in  this  program,  the 
APSI  B/Al  blades  will  be  tested  in  a  slot  angled  to  provide  normal  running  con¬ 
ditions.  This  will  result  in  a  higher  than  normal  incidence  angle  between  the 
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blade  and  the  injected  projectile  on  impact.  A  reduced  blade  rotational  speed 
will  be  employed  to  make  the  impact  normal  energy  transfer  more  nearly  the 
same  as  for  an  impact  at  takeoff  conditions.  Some  further  reduction  in  rota¬ 
tional  speed  may  be  necessary  to  permit  operation  with  the  existing  injection 
mechanism;  it  has  not  been  used  for  rotational  speeds  above  10,000  rpm. 

The  F101  blades  will  also  be  tested,  in  an  unrestrained  state  at  the  tip, 
in  a  disk  with  slots  angled  for  normal  running  conditions.  This  again  will 
result  in  a  higher  than  normal  incidence  angle  between  the  blade  and  the  in¬ 
jected  projectile  on  impact.  As  with  the  APS1  blades,  the  blade  rotational 
speed  will  be  reduced  to  provide  an  impact  normal  energy  transfer  that  would 
be  the  same  as  for  an  impact  at  takeoff  condition. 

No  problems  will  result  from  these  different  rotational  speeds  and  blade 
angles  since  the  damage  predictions  of  Tasks  V  and  VIII  will  be  made  using 
the  same  speeds  and  angles  as  for  these  test  conditions.  The  conditions  shown 
are  for  simulated  takeoff  conditions.  Blade  impact  locations  are  30  and  70% 
span.  The  size  of  the  bird  or  ice  slab  for  the  first  impact  at  each  location 
will  be  analytically  selected  to  be  just  below  the  blade-damage  condition. 

The  projectile  size  for  the  second  impact  at  a  location  will  be  based  on  ana¬ 
lytical  predictions  for  a  level  slightly  above  the  threshold  for  local  damage. 
Three  blades  of  each  configuration  will  be  used  for  the  bird/30%  span,  bird/70% 
span,  and  ice  slab/70%  span  tests.  The  fourth  blade  of  each  configuration  will 
be  held  in  reserve.  The  impact  conditions  selected  for  the  testing  of  these 
fourth  blades  will  be  determined  by  the  results  of  the  first  six  tests  of  each 
blade  configuration.  All  impacts  will  be  recorded  by  three  high-speed  cameras 
mounted  at  three  different  angles  to  obtain  detailed  visualization  of  the  im¬ 
pact  sequence  and  to  determine  the  resulting  blade  deflections. 

Each  test  blade  will  be  subjected  to  a  thorough  failure  analysis.  The 
failure  analysis  will  consist  of  visual  inspection  documented  with  photographs 
and  supplemented  by  dimensional  and  dye-penetration  or  Magnaflux  inspection. 
Also,  frequencies  in  the  first  four  modes  will  be  measured  for  any  blades  with 
local  impact  damage.  To  assess  the  significance  of  ther  damage,  selected  dam¬ 
aged  blades  will  be  subjected  to  bench  fatigue  tests  in  the  first-flexural  or 
first-torsional  mode  and  compared  to  available  or  control  fatigue-test  data 
for  these  blades. 


11.3  IMPACT  TESTING 


11.3.1  J79  Compressor  Stage  1  Blade  Testing 

The  J79  compressor  Stage  1  steel  blades  have  been  single-blade  tested,  in 
the  GE  Evendale  Cell  63  whirligig  facility,  using  the  disk  and  setup  shown  in 
Figure  60.  The  facility  is  comprised  of  an  electric  motor,  clutch,  protective 
shrouding,  bell  jar,  and  injection  mechanism.  The  J79  test  disk  has  a  number 
of  slots  with  different  setup  angles  so  that  the  blade  impact  angle  can  be  se¬ 
lected  to  suit  the  desired  impact  conditions.  The  bird-to-b lade  incidence 
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angle  was  set  to  33°  for  70%  radial-span  location  tests  and  47°  for  30%  radial- 
span  location  tests.  Rotor  speed  was  8243  rpm  for  the  70%  span  tests  and  8904 
rpm  for  the  30%  span  tests  and  provided  bird-to-blade  relative  velocities  of 
809  and  580  feet  per  second.  These  angles  and  speeds  simulate  an  injection  at 
takeoff  condition  and  are  the  same  as  would  occur  if  the  bird  were  shot  at  the 
engine  with  a  velocity  of  250  feet  per  second.  The  bird  was  simulated  by  a 
cylinder  made  of  a  microballoon  gelatin  mix. 

Figure  62  shows  the  J79  blade  and  disk  setup  and  the  injector  mechanism 
used  to  project  the  bird  into  the  path  of  the  blade.  Good  strain  data  and 
high-speed-movie  data  are  available  for  these  tests  and  are  currently  being 
ana  Ivzed . 


11.3.2  APSI  Fan  B/Al  Blade  Testing 

The  APSI  fan  B/Al  blades  will  be  single-blade  tested  in  the  GE  Evendale 
Cell  63  whirligig  facility  using  the  same  disk  and  setup  as  used  for  the  J7S 
blade.  The  APSI  blades,  however,  will  be  set  up  in  the  normal-blade-running 
angular  position  because  the  existing  test  disk  is  set  up  for  this  angle  for 
APSI  blade  testing.  The  plan  is  to  test  the  blade  at  a  rotor  speed  of  11,521 
rpm  for  the  70%  span  tests  and  at  13,846  rpm  for  the  30%  span  tests.  This 
will  provide  a  bird/blade  incidence  angle  of  27°  and  bird-to-blade  relative 
velocity  of  922  feet  per  second  for  the  70%  radial-span  location  testing  and 
incidence  angle  of  49°  and  bird/blade  relative  velocity  of  815  feet  per  seccnd 
for  30%  radial-span  location  testing.  These  conditions  provide  higher  inci¬ 
dence  angle  and  lower  relative  velocities  than  would  be  expected  for  normal 
takeoff  conditions  of  this  high-tip-speed  blade.  These  test  conditions  appear 
reasonable,  however,  for  the  following  reasons: 

1.  The  resulting  pro ject ile-to-blade  normal  energy  transfer  would  be 
approximately  the  same  as  for  an  impact  involving  a  lower  incidence 
angle  and  a  higher  relative  velocity. 

2.  Damage  predictions  will  be  made  for  the  same  relative  velocity  and 
incidence  angle  as  used  for  the  testing.  Predicted  bird  size  for 
no  damage,  threshold-minor  local  damage,  and  threshold-major  damage 
will  be  made. 


11.3.3  F101  Fan  Stage  1  Blade  Testing 

The  F101  fan  Stage  1  tip-shrouded  titanium  blades  will  be  single-blade 
tested  with  the  tip  shrouds  remaining  unsupported.  The  setup  will  be  similar 
to  that  used  for  the  J79  and  APSI  blades  except  that  the  F101  engine  disk  will 
be  set  up  with  a  soft-mount  arrangement  in  GE  Evendale  Cell  62  (Figure  61). 
Since  the  engine  disk  has  all  slots  arranged  for  normal  running,  the  bird-to- 
blade-impact  incidence  angle  will  be  higher  for  single-blade  testing  because 
of  the  zero  axial  speed  of  the  projectile. 
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The  rotational  speeds  will  be  taken  as  5t>42  rpm  for  7U%  span  tests  and 
6317  rpm  for  30%  span  tests  to  provide  the  same  normal  energy  for  impact  as 
a  bird  ingested  at  takeoff  conditions  while  accounting  for  the  differences 
in  incidence  angles.  A  projectile  velocity  of  200  feet  per  second  and  a  fan 
rotational  speed  of  7500  rpm  is  considered  to  be  a  typical  takeoff  condition 
for  this  engine  and  was  used  to  determine  the  above  test  rotational  speeds. 

H .4  WHIRLIGIG  FACILITY 

The  test-rig  facilities  (whirligig)  to  be  used  for  this  effort  were  de¬ 
veloped  to  closely  simulate  impact  conditions  that  occur  in  turbofan  engines. 

A  rig  typically  consists  of  a  horizontal  spindle  shaft,  electric  motor,  gear¬ 
box,  shrouding,  and  a  bell  jar  (Figure  63).  The  electric  motor  is  a  constant- 
speed  unit  with  a  dynamic,  variable-speed  output  and  magnetic  clutch.  Power 
is  transmitted  through  a  gearbox  and  a  horizontal  spind’e  shaft  to  the  rotor. 
Shrouding  is  provided  around  the  rotor  to  reduce  air  turbulence  and,  thus, 
heat  generation  and  power  requirements.  An  environmental  chamber  (belljar) 
is  soft-mounted  to  the  facility  shroud  to  provide  a  helium  atmosphere  to  re¬ 
duce  the  power  requirement  and  temperatures,  to  provide  a  housing  for  the  in¬ 
jection  mechanism,  and  to  enable  high-speed  movies  to  be  taken  of  the  impact 
sequence . 

The  soft-mounted  vehicle  setup  to  be  used  for  the  F101  blade  testing  is 
shown  in  Figure  61.  It  is  basically  a  standard  TF34  engine  fan  package  with 
some  additions  and  changes.  The  entire  vehicle  is  soft-mounted  and  is  sup¬ 
ported  by  spring-type  hangers  and  dampers.  The  test  disk  is  an  F101  Stage  1 
rotor  disk  adapted  to  the  setup  through  a  special  drive  cone.  The  rotor  is 
driven  by  a  flexible  coupling  and  slave  adapting  shaft.  The  extremely  high 
rotor  vibrations  resulting  from  loss  of  significant  portions  of  a  blade  dur¬ 
ing  impact  testing  are  reduced  by  the  spring-and-damper  mounting.  The  flexi¬ 
ble  drive  coupling  prevents  damage  to  the  facility  drive  spindle  and  bearings 
that  might  otherwise  result  from  large,  fan-rotor  displacements. 
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12.0  TASK  X  -  FULL-STAGE  RESPONSE  PREDICTION 


12.1  INTRODUCTION 


The  purpose  of  this  task  is  to  predict  the  initial  response  of  a  full 
set  of  blades  by  applying  the  final  analytical  models  of  Task  II,  the  final 
loading  models  of  Task  III,  and  the  final  damage  models  of  Task  VIII. 


12.2  ANALYSIS  APPROACH 

The  principal  difference  between  this  analysis  and  those  made  in  Task  V 
for  analysis  of  the  single-blade  rotating  tests  of  Task  IX  is  the  interactive 
influence  of  multiple  airfoils.  The  major  interactive  mechanism  occurs  when 
the  airfoils  are  interconnected  with  part-span  shrouds  or  tip  shrouds. 

The  effect  of  the  concentrated  stiffness  represented  by  the  shroud  will 
be  included  in  the  first-  and  second-level  analyses  of  Task  II  through  the 
use  of  spring  elements.  This  modeling  approach  will  be  modified  and  refined 
as  necessary  based  on  the  analyses  and  test  results  of  Tasks  V,  VI,  and  IX. 
The  results  of  that  effort  will  be  utilized  in  Task  X  to  make  the  full-stage 
response  prediction. 

An  additional  interaction  between  airfoils  can  take  place  through  the 
dovetails  and  disk,  but  this  is  greatly  dependent  upon  the  relative  stiffness 
of  the  disk  rims,  the  solidity  of  the  stage,  the  size  of  the  foreign  object, 
and  the  proximity  of  the  impact  to  the  root  of  the  airfoil.  In  general,  the 
interaction  at  the  dovetail  is  usually  much  less  than  the  interaction  at  the 
shrouds . 

The  cost  for  running  the  second-level  model  (COMET)  for  far-field  analy¬ 
sis  would  be  prohibitive  because  of  the  small  time  step  required.  We  will 
use  the  first- level  model  (NOSAPM)  for  far-field  analysis  through  the  restart 
capability  of  this  program,  and  both  C0MET-F0D  and  NOSAPM  will  be  used  for 
local- damage/ smal 1- time- per iod  anal ysi s . 


12.3  Response  Predictions 

Predictions  will  be  made  for  the  F101  shrouded  titanium  blades.  The  pre¬ 
dictions  will  span  the  range  of  expected  foreign  objects  (small  birds,  large 
birds,  and  ice)  and  will  include  relative  severities  ranging  from  minor 
damage  to  the  most  severe  engine  and  aircraft  ingestion  condition,  as  identi¬ 
fied  in  Task  I,  which  can  be  expected  to  result  in  severe  airfoil  damage. 

These  predictions  will  be  used  to  select  the  impact  parameters  for  tne  Task 
XI  full-scale  tests;  the  results  of  these  tests  will  be  used,  in  turn,  to 
make  the  final  correlative  adjustments  in  the  models  as  may  be  required. 
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13.0  TASK  XI  -  FULL-SCALE  IMPACT  TESTS 


13.1  INTRODUCTION 

The  object  of  Task  XI,  "Full-Scale  Impact  Tests,"  is  to  validate  the 
Task  X  prediction  of  how  a  full  stage  of  shrouded  blades  will  response  to  im¬ 
pact  loadings.  The  validation  will  be  accomplished  by  impact  testing  a  full 
stage  of  F101  tip-shrouded  blades  in  an  existing  whirligig  facility. 

The  tests  will  be  conducted  under  test  conditions  that  duplicate  the 
analytical  impact-parameter  values  used  for  prediction  purposes  in  Task  X. 
Enough  blades  will  be  instrumented  to  provide  the  response  data  needed  for 
the  validation.  The  use  of  duplicate  test  and  analytical  impact  parameters 
will  permit  a  direct  correlation  of  the  analytical  and  experimental  results. 

Activities  conducted  under  this  task  will  be  carefully  coordinated  with 
the  experimental  and  analytical  studies  of  Tasks  V,  VI,  IX,  and  X. 


13.2  TEST  PLAN 


Impact  tests  will  be  conducted  on  a  full  stage  of  F101  fan  blades  in  an 
existing  whirligig  facility.  The  F101  fan  Stage  1  tip-shrouded  titanium  blade 
was  selected  in  Task  V  for  single-blade  impact  testing  in  Task  IX  and  the 
full-stage  testing  of  this  Task  IX.  The  blades  will  be  tested  using  an  FlUl 
engine  disk  setup  in  a  soft-mount  TF34  vehicle  in  General  Electric  Evendale 
Cell  No.  62.  Figure  61  shows  the  test  setup. 

Two  tests  will  be  conducted.  In  the  first  test,  the  blades  will  be  im¬ 
pacted  by  a  real  bird;  in  the  second  test,  they  will  be  impacted  by  an  ice 
slab.  A  summary  of  the  impact-test  conditions  is  presented  in  Table  16.  Ro¬ 
tational  speed  and  projectile  velocities  will  be  set  to  simulate  ingestions 
at  takeoff  conditions.  For  both  tests,  the  rotational  speed  will  be  7300  rpm, 
and  the  impacts  will  be  at  70%  span.  The  bird  will  be  shot  into  the  blades 
at  a  velocity  of  200  ft/sec  and  the  ice  slab  at  75  ft/sec.  For  each  test, 
nine  blades  will  be  instrumented  with  three  gages  apiece.  This  is  sufficient 
instrumentation  to  record  the  responses  of  the  impacted  and  adjacent  blades. 

The  projectile  size  for  each  of  the  tests  will  be  determined  in  Task  X 
and  will  be  based  on  an  analytically  predicted  low  level  of  blade  damage  which 
slightly  exceeds  the  threshold  for  local  damage.  Both  impact  tests  will  be 
recorded  by  three  high-speed  cameras  mounted  from  three  different  angles  to 
obtain  detailed  views  of  the  impact  sequence  and  to  determine  the  resulting 
blade  deflections. 

All  damaged  blades  will  be  subjected  to  a  thorough  failure  analysis.  The 
failure  analysis  will  consist  of  visual  inspection  documented  with  photographs 
and  supplemented  by  dimensional  and  magnetic-particle  inspection. 


117 


13.3  IMPACT  TESTING 


13.3.1  F101  Fan  Stage  I  Blades  -  Bird-Ingestion  Test 

A  full  stage  of  F101  fan  Stage  1  tip-shrouded  titanium  blades  will  be 
subjected  to  a  single  bird  impact.  The  setup  will  be  similar  to  the  setup  foi 
the  single-blade  testing  in  Task  IX,  the  major  difference  being  that  the  disk 
will  contain  a  full  set  of  50  blades.  The  F101  disk  will  be  set  up  with  a 
soft-mount  arrangement  in  General  Electric  Evendale  Cell  63  (Figure  61). 

Since  an  engine  disk  will  be  used,  and  the  projectile  will  be  shot  into  the 
rotor  at  the  same  speed  expected  for  the  takeoff  condition  being  simulated, 
the  rotor  speed  will  also  be  the  same  as  for  actual  engine  operation.  This 
will  provide  conditions  identical  to  those  expected  for  a  similar  bird  inges¬ 
tion  during  engine  operations.  The  test  conditions  are  defined  in  Table  16. 


13.3.2  F101  Fan  Stage  1  Blades  -  Ice-Slab-Ingestion  Test 

The  ice-s Lab-inges t ion  test  will  follow  the  previously  described  bird- 
ingestion  test.  A  full  stage  of  F101  blades  will  be  tested.  The  more  severly 
damaged  blades  from  the  first  test  will  be  replaced  before  this  test  is  con¬ 
ducted.  The  point  of  impact  for  this  test  will  again  be  at  the  blade  70%  span 
location  but  will  be  targeted  to  occur  at  a  location  180°  circumferentially  on 
the  rotor  from  the  first  test  impact.  The  same  basic  test  setup  will  be  em¬ 
ployed,  and  test  conditions  will  be  selected  to  simulate  an  ice-slab  ingestion 
at  takeoff  conditions  (Table  16).  The  ice  slab  will  be  an  ice  spear  of  approx¬ 
imately  2  inch  by  2  inch  by  2  inch  triangular  cross  section  with  a  length 
compatible  with  the  mass  predicted  to  cause  an  appropriate  level  of  damage  in 
the  blades. 


13.4  WHIRLIGIG  FACILITY 


The  test-rig  facilities  (whirligig)  to  be  used  for  this  effort  were  de¬ 
veloped  to  closely  simulate  impact  conditions  that  occur  in  turbofan  engines. 

A  rig  typically  consists  of  a  horizontal  spindle  shaft,  electric  motor,  gear¬ 
box,  shrouding,  and  a  bell  jar  (Figure  63).  The  electric  motor  is  a  constant- 
speed  unit  with  a  dynamic,  variable-speed  output  and  magnetic  clutch.  Power 
is  transmitted  through  a  gearbox  and  a  horizontal  spindle  shaft  to  the  rotor. 
Shrouding  is  provided  around  the  rotor  to  reduce  air  turbulence  and,  thus, 
reduce  heat  generation  and  power  requirements.  An  environmental  chamber  (bell 
jar)  is  soft-mounted  to  the  facility  shroud  to  provide  a  helium  atmosphere  to 
reduce  the  power  requirement  and  temperatures,  and  to  enable  high-speed  movies 
to  be  taken  of  the  impact  sequence. 

The  vehicle  setup  to  be  used  for  the  F101  blade  testing  is  shown  in  Fig¬ 
ure  61.  It  is  basically  a  standard  TF34  engine  fan  package  with  sonic  additions 
and  changes.  The  entire  vehicle  is  soft-mounted  and  is  supported  by  spring- 
type  hangers  and  dampers.  The  test  disk  is  an  F101  Stage  1  rotor  disk  adapted 
to  the  setup  through  a  special  drive  cone.  The  rotor  is  driven  by  a  flexible 
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Table  16.  Planned  Test  Conditions  for  Full-Stage  Whirligig  Impact  Tests 


•  Simulated  Takeoff 

•  F101  Stage  1  Fan 


Test  No.  1 

Test  No .  2 

Blade  Type 

Tip-Shrouded 

Tip- Shrouded 

Blade  Material 

Titanium  8-1-1 

Titanium  8-1-1 

Number  of  Blades 

50 

50 

Number  of  Blades 
Instrumented 

9 

9 

Projectile  Type 

Real  Bird 

Ice  Slab 

Projectile  Speed 

200  ft/sec 

75  ft/sec 

Rotor  Speed 

7500  rptn 

7500  rpm 

Blade  Tip  Speed 

1468  ft/sec 

1468  ft/sec 

Impact  Location 

70%  Span 

70%  Span 

Incidence  Angle 

25° 

31° 

Relative  Velocity 

1263  ft/sec 

1245  ft/sec 

Basis  for  Projectile 
Size 

Threshold  Damage 

Threshold  Damage 

coupling  and  sLave  adapting  shaft.  The  extremely  high  rotor  vibrations  result¬ 
ing  from  loss  of  significant  portions  of  a  blade  during  the  impact  test  are 
reduced  by  the  spr ing-and-damper  mounting.  The  flexible  drive  coupling  pre¬ 
vents  damage  to  the  facility  drive  spindle  and  bearings  that  might  otherwise 
result  from  lrrge,  fan-rotor  displacements. 

13.5  GAS  GUN  FACILITY 


The  gas  gun  mechanism  basically  consists  of  a  pressure  vessel,  a  rupture 
disk  assembly,  and  the  gun  barrel.  The  bird  is  contained  in  a  cylindrical 
carrier  that  is  propelled  down  the  gun  barrel  when  the  pressurized  gas  (helium) 
is  released.  The  carrier  is  retained  at  the  end  of  the  barrel,  and  the  bird 
is  shot  out. 

A  mylar  membrane  acts  as  the  diaphragm  of  the  rupture  disk  assembly. 
Rupture  of  the  diaphragm  occurs  when  the  mylar  is  weakened  by  heat.  A  typical 
gas  gun  shot-triggering  sequence  might  be  as  follows. 

After  the  desired  rotor  speed  is  reached,  a  high-speed  camera  is  acti¬ 
vated  manually.  A  relay  in  the  camera  trips  when  the  film  is  at  speed.  The 
relay  activates  the  flash-lighting  system  and  the  whirligig  triggering  system. 
After  a  preset  time  delay,  the  whirligig  triggering  system  activates  a  hot 
wire  that  melts  the  mylar  diaphragm,  thus,  "firing"  the  projectile. 

The  technique  of  synchronizing  the  bird  with  the  rotor  to  impact  a  par¬ 
ticular  segment  of  blades  is  accomplished  by  using  the  whirligig  triggering 
system  and  the  hot-wire  power-supply  unit.  The  hot-wire  power-supply  unit 
was  designed  and  built  specifically  for  the  purpose  of  rupturing  the  diaphragm 
rapidly  and  predictably.  The  whirligig  triggering  system  has  a  very  precise 
timing  circuit  that  can  be  varied  and  is  dialed  in  to  allow  for  the  operating 
time  of  the  gas  gun  system  and  to  account  for  the  angular  position  of  the 
blade  segment  relative  to  a  1/rev  indicator  on  the  rotor. 
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